10/07 03 THU 15:48 FAX 972 3 6127575 



G. E, EHRLICU 



g]018 



APPLICATION FOR PATENT 



Inventors: 



YoavKIMCIIY,RoniAMRAMI, Yona BOUSKILA, 
Udi ANTEBI, Nick SIDORENKO, Gal BEN-DAVID 
and Yoel Zilberstein 



5 



Title: 



RADIOACTIVE EMISSION DETECTOR EQUIPPED 
WITH A POSITION TRACKING SYSTEM 



10 



This is a continuation-in-part of U.S. Patent Application No. 1 0/343,792, filed 
on My 1 1, 200U which is a National Phase application of PCT/IL0I/00638» filed on 
July U, 2001, which is a continuation-in-part of U.S. Patent Application No. 

15 09/727,464 filed December 4, 2000, which is a continuation-in-part of U.S. Patent 
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FIELD AND BACKGROUND OF THE INVENTION 

The present invention relates to a radioactive emission probe equipped with a 

25 position tracking system. More particularly, the present invention relates to the 
functional integration of a radioactive emission probe equipped with a position 
tracking system as above with medical imaging modalities and (or) with guided 
minimally-invasive surgical instruments. The present invention is therefore useful for 
calculating the position of a concentrated radiophannaceutical in the body in positional 

30 context of imaged portions of the body, which information can be used, for example, 
for performing an efiScient minimally invasive sutgical procedure. The present 
invention further relates to a surgical instrument equipped with a position tracking 
system and a radioactive emission probe for fine in situ localization during resection 
and (or) biopsy procedures, which surgical instrument is operated in concert with other 

35 aspects of die invention. 

The use of minimally invasive surgical techniques has dramatically affected 
the methods and outcomes of surgical procedures. Physically cutting through tissue 
and organs to visually expose surgical sites m conventional "open siu'gical" 
procedures causes tremendous blunt tmuma and blood loss. Exposure of internal 

40 tissues and organs in this manner also increases the risk of infection. Traiuna, blood 
loss, and infection all combine to extend recovery times, increase the rate of 
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complications, and require a more intensive care and monitoring regiment. The result 
of such open surgical procedures is more pain and suffering, higher procedural cx)sts^ 
and greater risk of adverse outcomes. 

hi sharp contrast, minimally invasive surgical procedures cause little blunt 

S trauma or blood loss and minimize the risk of infection by maintaining the body's 
natural barriers to infection substantially intact. Minimally invasive surgical 
procedures result in faster recovery and cause fewer complications than conventional, 
open, surgical procedures. Minimally invasive surgical procedures, such as 
laparoscopic, endoscopic, or cystoscopic surgeries, have replaced more invasive 

10 surgical procedures in all areas of surgical medicine. Due to technological 
advancements in areas such as fiber optics, micro-tool fabrication, imaging and 
material science, the physician performing the operation has easier-to-operate and 
more cost effective tools for usie in minimally invasive procedures. However, there 
still exist a host of technical hurdles that hmit the efficacy and increase the difficulty 

15 of minimally invasive procedures, some of which were overcome by the development 
of sophisticated imaging techniques. As is fiirther detailed below the present invention 
offers a yet fiirther advantage in this respect. 

Radionuclide imaging is one of the most important applications of radioactivity 
in medicine. The purpose of radionuclide imaging is to obtain a distribution image of 

20 a radioactively labeled substance, e.g., a radiopharmaceutical, within the body 
following administraiion thereof to a patient. Examples of radiopharmaceuticals 
include monoclonal antibodies or other agents, e.g., fibrinogen or fluorodeoxyglucose, 

tagged with a radioactive isotope, e.g., ^^^technetium, ^^gallium, thallium, 

111 123 125 18 

indium, iodine, iodine and fluorine, which may be administered orally or 

25 intravenously. The radiophaimaceuticals are designed to concentrate in the area of a 
tumor, and the uptake of such radiopharmaceuticals in the active part of a timior, or 
other pathologies such as an inflaimnation, is higher and more rapid than in the tissue 
that neighbors the tumor. Thereafter, a radiation emission detector, typically an 
invasive detector or a gamma camera (see below), is employed for locating the 

30 position of the active area. Another application is the detection of blood clots with 
radiopharmaceuticals such as ACUTBCT from Nycomed Amersham for the detection 
of newly formed thrombosis in veins, or clots in arteries of the heart or brain, in an 
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emergency ox operating room. Yet other applications include radioimagiog of 
myocaxdial infarct using agents such as radioactive anti-rayosin antibodies^ 
radioimaging specific cell types using radioactively tagged molecules (also known as 
molecular imaging), etc. 

5 The distribution image of the radiophannaceuticai in and around a tumor, or 

another body stTuctiune> is obtained by recording the radioactive emission of the 
radiopharmaceutical with an external radiation detector placed at different locations 
outside the patient The usual preferred emission for such applications is that of 
gamma rays, which emission is in the energy range of ^proximately 20-511 KeV. 

10 When the probe is placed in contact with the tissue, beta radiation and positrons may 
also be detected. 

The first attempts at radionuclide "imaging" were in the late 1940*5. An anray 
of radiation detectors was positioned mechanically on a matrix of measuring points 
around the head of a patient. Alternatively, a single detector was positioned 

IS mechanically for separate measurements at each point on the matrix. 

A significant advance occurred in the early 1950's with the introduction of the 
rectihnear scanner by Ben Cassen. With this instrument, the detector was scanned 
mechanically in a predetennined pattern over the ar^ of interest. 

The first gamma camera capable of recording all points of the image at one 

20 time was described by Hal Anger in 1953. Anger used a detector comprised of a 
Nal(Tl) screen and a sheet of X-ray fitoi. In the late 1950's, Anger replaced the film 
screen with a photomuitipher tube assembly. The Anger camera is described in Hal 0. 
Anger, "Radioisotope camera in Hine GJ", Instrumentation in Nuclear Medicine, New 
York, Academic Press 1967, chapter 19. U.S. Patent No. 2,776,377 to Anger, issued 

25 in 1 957, also describes such a radiation detector ass«nbly. 

U.S. Patent No. 4,959,547 to Carroll et al. describes a probe used to map or 
provide imaging of radiation within a patient. The probe comprises a radiation 
detector and an adjustment mechanism for adjusting the solid angle through which 
radiation may pass to the detector, the solid angle being continuously variable. The 

30 probe is constructed so that the only radiation reaching the detector is that which is 
within the solid angle. By adjusting the solid angle from a maximum to a minimum 
while moving the probe adjacent the source of radiation and sensing the delected 
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radiation, one is able to locate the probe al the source of radiation. The probe can be 
used to determine the location of the radioactivity and to provide a point-by-point 
image of the radiation source or data for mapping the same. 

U.S. Patent No. 5;246,005 to Carroll et al. describes a radiation detector or 

« 

5 probe^ which uses statistically valid signals to detect radiation signals from tissue. The 
output of a radiation detector is a series of pulses, which are counted for a 
predeteimined amount of time. At least two count ranges are defmed by circuitry in 
the apparatus and the count range which includes die input count is determined. For 
each count range» an audible signal is produced which is audibly distrainable from the 

10 audible signal produced for every other count range. The mean values of each count 
range are chosen to be statistically different, e.g.» 1, 2, or 3 standard deviations, from 
the mean of adjacent lower or higher count ranges. The parameters of the audible 
signal, such as frequency, voice, repetition rate, and (or) intensity are changed for each 
count range to ptovide a signal which is discriminable firom the signals of any other 

IS count range. 

U.S. Pat^t No. 5,475,219 to Olson describes a system for detecting photon 
emissions wherein a detector serves to derive electrical parameter signals having 
amplitudes corresponding with the detected energy of the photon emissions and other 
signal generating events. Two comparator networks employed within an ^ergy 

20 window, which define a function to develop an ouq>ut, L, when an event-based signal 
ampUtude is equal to or above a threshold value, and to develop an output, H, when 
such signal amplitude additionally extends above an upper limit. Improved rehability 
and accuracy is achieved with a discriminator circuit which, in response to these 
ou^uts L and H, derives an event output upon the occurrence of an output L in the 

25 absence of an output H. This discriminator circuit is an asynchronous, sequential, 
fundamental mode discriminator circuit with three stable states. 

U.S. Patent Nos. 5,694,933 and 6.135,955 to Madden et al. describe a system 
and method for diagnostic testing of a structure within a patient's body that has been 
pnDvided with a radioactive imaging agent, e.g., a radiotracer, to cause the structure to 

30 produce ganuna rays, associated characteristic x rays, and a continuum of 
Compton-scatt^ed photons. The system includes a radiation receiving device, e.g., a 
hand-held probe or camera, an associated signal processor, and an analyzer. The 
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radiation receiving device is arranged to be located adjacent the body and the structure 
for receiving gamma rays and characteristic X-rays emitted from the structure and for 
providing a processed electrical signal representative thereof. The processed electrical 
signal includes a first portion representing the characteristic X-rays received and a 
5 second portion representing the gamma rays received. The signal processor removes 
flie signal corresponding to the Conqston-scattered photons from the electrical signal in 
the region of the ftill-energy gamma ray and the characteristic X-ray. The analyzer is 
arranged to selectively use the X-ray portion of the processed signal to provide 
near-field information about the structure^ to selectively use both the X-ray and the 

10 ganuna-ray portions of the processed signal to provide near*field and far-field 
information about the structure, and to selectively use the gamma-ray portion of the 
processed signal to provide extended field information about the structure. 

U.S. Patent No. 5,732,704 to Thurston et al, describes a method for identifying 
a sentinel lymph node located within a grouping of regional nodes at a lymph drainage 

15 basin associated with neoplastic tissue wherein a radiopharmaceutical is injected at the 
situs of the neoplastic tissue. This radiopharmaceutical migrates along a lymph duct 
towards the drainage basin containing the sentinel node. A hand-held probe with a 
forwardly disposed radiation detector crystal is maneuvered along the duct while the 
clinician observes a graphical readout of count rate amplitudes to determine when the 

20 probe is aligned with the duct The region containing the sentinel node is identified 
when the count rate at the probe substantially increases. Following surgical incision, 
the probe is maneuvered utilizing a sound output in connection with actuation of the 
probe to establish increasing count rate thresholds followed by incremental movements 
until the threshold is not reached and no sound cue is given to the surgeon. At this 

25 point of the maneuvering of the probe, the probe detector will be in adjacency with the 
sentinel node, which then may be removed. 

U.S. Patent No. 5,857,463 to Thurston et al. describes further apparatus for 
tracking a radiopharmaceutical present within the lymph duct and for locating the 
sentinel node within which the radiopharmaceutical has concentrated. A smaller, 

30 straight, hand-held probe is employed carrying two hand actuable switches. For 
tracking procedures, the probe is moved in an undulatory manner, wherein the location 
of the radiopharmaceutical^containing duct is determined by observing a graphics 
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readout When the region of the sentinel node is approached, a switch on the probe 
device is actuated by the surgeon to carry out a sequence of squelching operations until 
a small node locating region is defined. 

U.S. Palait Nos. 5^16,167 to Kramar et al. and 5,987,350 to Thurston describe 
S surgical probes wherein a heat-sieriiizahle and reusable detector component is 
combined with a disposable handle and cable assembly. The reusable detector 
component incoiporates a detector crystal and associated mountings along with 
preamplifier components. 

U.S. Patent No. 5,928,150 to Call describes a system for detecting emissions 

10 from a radiopharmaceutical injected within a lymph duct wherein a hand-held probe is 
utilized. When employed to locate sentinel lymph nodes, supplementary features are 
provided including a function for treating validated photon event pulses to detennine 
count rate level signals. The system includes a function for count-rate based ranging 
as well as an adjustable thresholding feature. A post-threshold amplification circuit 

15 develops Adl-scale aural and visual outputs. 

U.S. Patent Nos. 5,932,879 and 6,076,009 to Raylman et al. describe an 
intraoperative syst^ for preferentially detecting beta radiation over gamma radiation 
emitted from a radiopharmaceutical. The system has ion-implanted silicon 
charged-particle detectors for generating signals in response to received beta particles. 

20 A preamplifier is located in proximity to the detector filters and amplifies the signal. 
The probe is coupled to a processing unit for amplifying and filtering the signal. 

U.S. Patent No. 6,144,876 to Bouton describes a system for detecting and 
locating sources of radiation, with particular applicability to interoperative lymphatic 
mapping (EM) procedures. The scanning probe employed with the system performs 

25 with both an audible as well as a visual perceptive output. A desirable stability is 
achieved in the readouts from the system through a signal processing approach which 
establishes a floating or dynamic window analysis of validated photon event counts. 
This floating vrindow is defined between an upper edge and a lower edge. The values 
of these window edges vary during the analysis in response to compiled count sum 

30 values. In general, the upper and lower edges are spaced apart a value corresponding 
with about four standard deviations. 
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To compute these count sums, counts are collected over successive short scan 
intervals of SO milliseconds and the count segments resulting therefrom are located in 
a succession of bins within a circular buffer memory. The count sum is generated as 
the sum of the memory segment count values of a certain number of the bins or 
5 segm^ts of memory. Alteration of the floating window occurs when the count sum 
either exceeds its upper edge or falls below its lower edge. A reported mean, 
con4>uted with respect to the window edge that is crossed, is developed for each scan 
interval which, in turn, is utilized to derive a mean count rate signal. The resuhing 
perceptive ou^ut exhibits a desirable stability, particularly under conditions wherein 

10 the probe detector is in a direct confrontational geometry with a radiation source. 

U.S. Patent No. S,846,S13 teaches a system for detecting and destroying living 
tumor tissue within the body of a living being. The system is arranged to be used with 
a tumor localizing radiophannaceutical. The system includes a percutaneously 
insertable radiation detecting probe, an associated analyzer, and a percutaneously 

15 insertable tumor removing instrument, e.g., a resectoscope. The radiation detecting 
probe includes a needle unit having a radiation sensor component therein and a handle 
to which the needle unit is releasably mounted. The needle is arranged to be inserted 
through a small percutaneous portal into the patient's body and is movable to various 
positions within the suspected tumor to detect the pre.%nce of radiation indicative of 

20 cancerous tissue. The probe can then be removed and the tumor removing instrument 
inserted througji the portal to destroy and (or) remove the cancerous tissue. The 
instrument not only destroys the tagged tissue, but also removes it from the body of the 
being so that it can be assayed for radiation to confirm that the removed tissue is 
cancerous and not healthy tissue. A collimator may be used with the probe to establish 

25 the probe's field of view. 

The main limitation of the system is that once the body is penetrated, scanning 
capabilities are hmited to a translational movement along the line of penetration. 

An effective collimator for gamma radiation must be several mm in thickness 
and therefore an effective collimator for high energy gamma radiation cannot be 

30 engaged with a fine surgical instrument such as a surgical needle. On the other hand, 
beta radiation is absorbed mainly due to its chemical reactivity after passage of about 
0.2-3 mm through biological tissue. Thus, the system described in U.S. Patent No. 
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5,846,513 cannot efficiently employ high energy gamma detection because 
dixectionality will to a great extent be lost and it also cannot efficiently employ beta 
radiation because too high proximity to the radioactive source is required, whereas 
body tissue limits the degree of maneuvering the instrument. 

5 The manipulation of soft tissue organs requires visualization (imaging) 

techniques such as computerized tomography (CT), fluoroscopy pC-ray fluoroscopy), 
magnetic resonance imaging (MRI), optical endoscopy, mammography or ultrasound 
which distinguish the borders and shapes of soft tissue organs or masses. Over the c 
years, medical imaging has become a vital part in the early detection, diagnosis and 

10 treatment of cancer and other diseases. In some cases medical imaging is the first step 
in preventing the spread of cancer through early detection and in many cases medical 
imaging makes it possible to cure or eliminate the cancer altogether via subsequent 
treatment 

An evaluation of the presence or absence of tumor metastasis or invasion has 

15 been a major determinant for the achievement of an effective treatment for cancer 
patients. Studies have determined that about 30 % of patients with essentially newly 
diagnosed tumor will exhibit clinically detectable metastasis. Of the remaining 70 % 
of such patients who are deemed "clinically fi:^ee" of metastasis, about one-half are 
curable by local tumor therapy alone. However, some of these metastasis or even early 

20 stage primary tumors do not show with the imaging tools described above. Moreover 
often enough the most important part of a tumor to be removed for biopsy or surgically 
removed is the active, i.e., growing part, whereas using only conventional imaging 
caimot distinguish this specific part of a tumor from other parts thereof and (or) 
adjacent non affected tissue. 

25 A common practice in order to locate this active part is to mark it with 

radioactivity tagged materials generally known as radiopharmaceuticals, which arc 
administered orally or intravenously and which tend to concentrate in such areas, as 
the uptake of such radiopharmaceuticals in the active part of a tumor is higher and 
more rapid than in the neighboring tumor tissue. Thercafl^, a radiation emission 

30 detector, typically an invasive detector, is employed for locating the position of the 
active area. 
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Medical imaging is often used to build computer models which allow doctors 
to, for example, guide exact radiation in the treatment of cancer, and to design 
minimally-invasive or open surgical procedures. Moreover, imaging modalities are 
also used to guide surgeons to the target area inside the patient's body, in the operation 

5 room during the surgical procedure. Such procedures may include, for example, 
biopsies, inserting a localized radiation source for direct treatment of a cancerous 
lesion, known as brachytherapy (so as to prevent radiation damage to tissues near the 
lesion), injecting a chemotherapy agent into the cancerous site or removing a 
cancerous or other lesions, 

10 The aim of all such procedures is to pin-point the target area as precisely as 

possible in order to get the most precise biopsy results, preferably from the most active 
part of a timior, or to remove such a tumor in its entirety on the one hand with minimal 
damage to the surrounding, non affected tissues, on the other hand. 

However, in the current state of the prior ait this goal is yet to be achieved, as 

15 most of the common imaging modalities such as fluoroscopy, CT, MRI, 
mammognq)hy or ultrasound demonstrate the position and appearance of the entire 
lesion with anatomical modifications that the lesion causes to its surrounding tissue, 
without differentiating between the non-active mass from the physiologically active 
part thereof 

2U On the other hand, prior art radiation emission detectors and (or) biopsy 

probes, while being suitable for identifying the location of the radiation site, they leave 
something to be desired from the standpoint of facilitating the removal or other 
destruction of the detected cancerous tissue with minimum invasion of the patient. 

The combination of modalities, as is offered by the present invention, can 

25 reduce the margin of error in positioning such tumors. In addition, the possibility of 
demonstrating the position of the active part of a tumor superimposed on a scan from 
an imaging modality that shows the organ or tumor, coupled with the possibihty to 
follow a surgical tool in reference to the afflicted area during a surgical procedure will 
allow for a more precise and controlled surgical procedures to take place, minimizing 

30 the aforementioned problems. 

The present invention addresses these and other issues which are further 
elaborated hereinbelow, and offers the physicians and patients more reliable targeting, 
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that in turn will result in less invasive and less destructive surgical procedures and less 
cases of mistaken diagnosis. 

SUMMARY OF THE INVEPfTION 

4 

5 The present invention successfully addresses the shortcomings of the presently 

known configurations by providing a radioactive emission probe in communication 
with a position tracking system and the use thereof in a variety of systems and methods 
of medical imaging and procedures. Specifically, wide-{q>eture coUimation - 
deconvolution algorithms are provided, for obtaining a high-efficiency, high resolution 

10 image of a radioactivity emitting source, by scaiming the radioactivity emitting source 
with a probe of a wide-aperture collimator, and at the same time, monitoring the 
position of the radioactive emission probe, at very fine time intervals, to obtain the 
equivalence of fine-ap^ure coUimation. The blurring effect of the wide aperture is 
then corrected mathematically. Furthermore* an imaging method by depth calculations 

15 is provided, based on the attofiuation of photons of different energies, which are 
emitted firom the same source, coupled with position monitoring. 

The present invention has many other applications in the direction of 
Aerapeutics, such as, but not limited to, implanting brach3rtherapy seeds, ultrasound 
microwave radio-firequency cryotherapy and localized radiation ablations. 

20 Implementation of the methods and systems of the preset invention involves 

performing or completing selected tasks or steps manually, automatically, or a 
combination thereof. Moreover, according to actual instrumentation and equipment of 
preferred embodiments of the methods and systems of die present invention, several 
selected steps could be implemented by hardware or by software on any operating 

25 system of any firmware or a combination thereof. For example, as hardware, selected 
steps of the invention could be implemented as a chip a circuit. As software, selected 
stq)S of the invention could be implemented as a plurality of software instructions 
being executed by a computer using any suitable algorithms, hi any case, selected 
steps of the method and system of the invention could be described as being performed 

30 by a data processor, such as a computing platform for executing a plurality of 
instructions. 
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BRIEF DESCRimON OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to 
the accompanying drawings. With specific reference now to the drawings in detail, it 
is stressed that the particulars shown are by way of example and purposes of 

S illustrative discussion of the preferred embodiments of the present invention only, and 
are presented in the cause of providing what is believed to be the most useful and 
readily understood description of the principles and conceptual aspects of the 
inv^tion. In this regard, no attempt is made to show structural details of the invention 
in more detail than is necessary for a fundamental understanding of the invention, the 

10 description taken with the drawings making apparent to those skilled in the art how the 
several forms of the invention may be embodied in practice. 
In the drawings: 

FIG. 1 is a black box diagram of a system according to the teachings of the 
present invention; 

15 FIG. 2 is a perspective view of an articulated arm which serves as a position 

tracking system shown carrying a radioactive emission probe in accordance with the 
teachings of the present invention; 

FIG. 3 is a schematic depiction of a radioactive emission probe carrymg a pair 
of three coaxially aligned accelerometers which serve as a position tracking system in 
20 accordance with the teachings of the present invention; 

FIG. 4 is a schematic presentation of a radioactive emission probe 
communicating with yet another type of a position tracking system in accordance with 
the teachings of the present invention; 

FIG. 5 is a simplified cross-sectional view of a narrow or wide angle 
25 radioactive emission probe used to implement an embodiment of the present invention; 

FIG- 6 is a presentation of a scaiming protocol which can be effected with the 
detector of Figure 5; 

FIG. 7 is a simplified cross-sectional view of a spatially saisitivc radioactive 
emission probe, e.g., a gamma camera, used to implement another ^bodiment of the 
30 present invention; 

FIG. 8 is a presentation of a scanning protocol which can be effected with the 
detector of Figure 7; 
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FIG. 9 demonstrates a system in accordance with the teachings of the present 
invention which employs four position tracking systems for co-tracking the positions 
of a patient, a radioactive emission probe, an imaging modality and a surgical 
instrument; 

5 FIG. 10 demonstrates the use of a pair of radiation emission detectors 

connected therebetween via a connector, preferably a flexible connector or a flexible 

connection to the connector according to the present inv^tion; 

FIG. 1 1 is a schematic diagram of a surgical instrument and accompanying 

system elements according to the teachings of the present invention; 
10 FIG. 12 is a simplified pictorial illustration of an imaging system constructed 

and operative in accordance with a preferred embodiment of the present invention^ 

including a radiation probe and position sensor, position tracking system, medical 

imaging system and coordinate registration system; 

FIG. 13 is a simplified pictorial illustration of a single dimension image 
15 formation with a nuclear radiation probe attached to a position tracking system of the 

system of Figure 12, in accordance with a prefened embodiment of the present 

invention; 

FIG. 14 is a simplified pictorial plot of detecting a radiation point source with 
the nuclear radiation probe of the syst^ of Figure 12, without further processing, in 
20 accordance with a preferred embodiment of the present invention; 

FIG. 15 is a simplified flow chart of an averaging algorithm used in the 
imaging system of Figure 12, in accordance with a preferred embodiment of the 
present invention; 

FIG. 16 is a simplified pictorial plot of detecting a radiation point source with 
25 the nuclear radiation probe of the system of Figure 12, with averaging processing, in 
accordance with a preferred embodiment of ttie present mvention; 

FIGS. 17 and 18 are simplified pictorial illustrations of hot cross and hot bar 
phantom images, respectively, of images produced by a gamma radiation probe of the 
system of Figure 12; 

30 FIG. 19 is a simplified flow chart of a minimizing algorithm used in the 

imaging system of Figure 12, m accordance with a preferred embodiment of the 
present inventira; 
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FIG. 20 is a simplified pictorial plot of detecting a radiation point source with 
the nuclear radiation probe of the system of Figure 12, with minimizing processing, in 
accordance with a preferred embodiment of the present invention; 

FIG. 21 is a simplified pictorial illustration of an image reconstruction system 

« 

5 constructed and operative in accordance with a preferred embodiment of the present 
invention, which produces a combined image made up of medical images, the position 
of the peak radiation location and the location of a therapeutic instrument; 

FIG. 22 is a simphfied flow chart of a radiation map reconstruction algorithm, 
in accordance with a preferred embodiment of the present invention; 
10 FIGs. 23A and 23B are illustrations of radiolabeled patterns observed in 

images produced by the system of the invention and by a conventional gamma camera, 
respectively, of an autonomous adenoma of a thyroid; 

FIGs. 24A and 24B are illustrations of radiolabeled patterns observed in 
images produced by the system of the invention and by a conventional gamma camera, 
IS respectively, of suspected Paget ^s disease of a humerus; 

FIGs. 2SA and 25B are illustrations of radiolabeled patterns observed in 
images produced by the system of the invoition and by a conventional gamma camera, 
respectively, of chronic osteomyelitis; 

FIGs. 26A and 26B ai^ illustrations of radiolabeled patterns observed in 
20 images produced by the system of the invention and by a conventional ganmia camera, 
respectively, of skeletal metastasis from medutloblastoma; 

FIGs. 27A-27I demonstrate the operation of an algorithm provided by the 
present inv^tion for estimating the distribution of radiation sources in a control 
volume; 

25 FIGs. 28A - 28F schematically illustrate a handheld probe, in accordance with 

preferred embodim^ts of the present invention; 

FIGs. 29A - 29B schematically illustrates the manner of calibrating the 
handheld probe of FIGs 28A - 28F, in accordance with a preferred embodiment of die 
present invention; 

30 FIG. 30. schematically illustrates the manner of synchronizing event and 

position readings of the handheld probe of FIGs 28A - 28F, in accordance with a 
preferred embodiment of the present inv^tion; 
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FIGs. 31 A - 31C describe a spatial resolution bar-phantom Icsl of a 
radioactive emission probe, in accordance with a preferred embodiment of the present 
invention; 

FlGs. 32A - 32D describe a spatial resolution bar-phantom test of a prior-art 

5 probe; 

FIGs. 33A -33B illustrate the energy resolution of a single pixel of a 
radioactive emission probe» in accordance with the present invention; 

FIGs. 34A - 34C illustrate endoscopic radioactive emission probes; 
FIG. 35 illustrate a method of calculating the depth of a radiation source; and 
10 FIG. 36 illustrate a two-dimensional image of the radioactivity emitting 

source, produced by a free-hand scanning of a cancerours prostate gland, ex vivo. 

DESCRIPTION OF THE PREFERRED EMBOPIMENTS 

The present invention is of a radioactive emission probe in communication 

15 with a position tracking system and the use thereof in a variety of systems and methods 
of medical imaging and procedures. Specifically, wide-2^ture coUimation - 
deconvolution algorithms are provided, for obtaining a high-efiRciency, high resolution 
image of a radioactivity ooaitting source, by scanning the radioactivity emitting source 
with a probe of a wide-aperture collimator, and at the same time, monitoring the 

20 position of the radioactive emission probe, at very fine time intervals, to obtain the 
equivalence of fine-aperture collimation. The blurring effect of the wide aperture is 
then corrected mathematically. Furthermore, an imaging method by depth calculations 
is provided, based on the attentuation of photons of different energies, which are 
emitted from the same source, coupled with position monitoring. 

25 The principles and operation of the present invention may be better understood 

with reference to the drawings and accompanying descriptions. 

Before explaining at least one embodiment of the invention in detail, it is to be 
understood that the invention is not limited in its application to the details of 
construction and the arrangement of the components set forth in the following 

30 description or illustrated in the drawings. The invention is capable of other 
embodiments or of being practiced or earned out in various ways. Also, it is to be 
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understood that the phraseology and tenninology employed herein is for the purpose of 
description and should not be regarded as limiting. 

Functional imaging, or the use of radioactive materials to tag physiologically 
active tissue within the body of a patient, for detennining the tissue*s localization and 

■ 

5 d«narcation by radioactive emission probes has been disclosed in the medical 
literature for at least forty years. Functional imaging shows the metabolic activity of 
body tissue, since dying or damaged body tissue absorbs radiopharmaceuticals at a 
different rate from a healthy tissue. The functional image may be used for example, to 
study cardiac rhythm, or respiratory rhythm. However, a functional image may not 

10 show structural, or anatomic details. 

Significant developments in the locaUzation and demarcation of tissue bearing 
radioactive isotope tags for diagnostic and (or) therapeutic purposes have occurred 
since that time. In fact, it is now becoming an established practice in the diagnosis and 
(or) treatment of certain diseases, e.g., cancer, blood clots, myocardial infarct and 

IS abscesses, to introduce monoclonal antibodies or other agents, e.g., fibrinogen, 
fluorodeoxyglucose labeled with a radioactive isotope (e.g., ^^Technetium, 
^Gallium, ^^^ThalUum. ^^*Indium, '^^lodine, ^ ^Fluorine and ^^^lodine) into the body of 
the patient. Such radiopharmaceuticals tend to localize in particular tissue or celt type, 
whereas uptake or binding of the specific radiopharmaceutical is increased in more 

20 '"physiologically active" tissue such as tiie active core of a cancerous tissue, so that the 
radiation emitted following nuclear disintegrations of the isotope can be detected by a 
radiation detector to better allocate the active portion of a tumor. Such radiation may 
be, for example, a, |) , p and (or) y radiation. 

In another type of applications radioactive substances are used to determine the 

25 level of flow of blood in blood vessels and the level of perfusion thereof into a tissue, 
e.g., coronary flow and myocardial perfusion. 

Referring now to the drawings, Figure 1 illustrates a system for calculating a 
position of a radioactivity emitting source in a system-of-coordinates, in accordance 
with the teachings of the present invention, which system is referred to hereinbelow as 

30 system 20. 

System 20 includes a radioactivity emission detector 22. System 20 according 
to the present invention further includes a position tracking system 24. System 24 is 
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connected to and (or) communicating with radioactive cnussion probe 22 so as to 
monitor the position of detector 22 in a two- or three-dimensional space defined by a 
system-of-coordinates 28 in two, three or more, say four, five or preferably six 
degrees-of-lrcedom (X, Y, Z, p, 9 and cp). System 20 further includes a data processor 

5 26. Data processor 26 is designed and configured for receiving data inputs from 
position tiackiiig system 24 and JQrom radioactive emission probe 22 and» as is fiuther 
detailed below, for constructing the image of the radioactivity emitting souixre in 
system-of-coordinates 28. As shown in Figure 10, a pair (or more) of detectors 22 
connected therebetween via a physical connector, each of detectors 22 is position 

10 tracked, can be used for constructing the image of the radioactivity emitting source in 
system-of-coordinates 28. If more than a single detector 22 is used, detectors 22 are 
preferably connected there between via a connector 29. Connector 29 is preferably 
flexible. In the alternative, the connections of detectors 22 to connector 29 provide the 
required flexibility. 

15 It will be appreciated that system 20 of radioactivity emissim detector 22 and 

position tracking system 24 is inherently different from known SPECT and PET 
imaging systems, as well as from other imaging systems such as X-ray, 
Mammography, CT, and MRI, since the motion of detector 22 of the present invention 
is not limited to motion within a predetermined track or tracks, with a respect to an 

20 immovable gantry. Rather^ detector 22 of the present invention is adapted for a 
variable-course motion, which may be for example, free-hand scanmng, 
variable-course motion on a linkage system, motion within a body lumen, endoscopic 
motion through a trucar valve, or another form of variable-course motion. 

Position tracking systems per se are well known in the art and may use any one 

25 of a plurality of approaches for the determination of position in a two- or 
three-dimensional space as is defined by a system-of-coordinates in two, three and up 
to six degrees-of-frccdom. Some position tracking systems employ movable physical 
connections and appropriate ihovemeni monitoring devices (e.g., potentiometers) to 
keep track of positional changes. Thus, such systems, once zeroed, keep track of 

30 position changes to thereby determine actual positions at all times. One example for 
such a position backing system is an articulated arm. 
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Figure 2 shows an articulated arm 30 which includes six arm members 32 and a 
base 34, which can therefore provide positional data in six degrees-of-freedom. 
Monitoring positional changes may be efTected in any one of s vera! different ways. 
For example, providing each ann member 32 with, e.g., potentiometers or optical 

S encoders 38 used to monitor the angle between adjacent arm members 32, to thereby 
monitor the angular change of each such ann member with respect to adjacent arm 
members, so as to detemiine the position in space of radioactive emission probe 22, 
which is physically connected to articulated ann 30. 

As is shown in Figure 3 other position tracking systems can be attached 

10 directly to radioactive emission probe 22 in order to monitor its position in space. An 
example of such a position tracking system is an assortment of three triaxially (e.g., 
co-orthogonally) oriented accelerometers 36 which may be used to monitor the 
positional changes of radioactive emission probe 22 with respect to a space. A pair of 
such assortments, as is specifically shown in Figure 3, can be used to determine the 

15 position of detector 22 in six-degrees of freedom. 

As is shown in Figures 4 and 10, other position tracking systems re-determine a 
position iirespective of previous positions, to keep track of positional changes. Such 
systems typically employ an array of receivers/transmitters 40 which are spread in 
known positions in a system-of-coordinates and transmitt^(syreceiver(5) 42, 

20 respectively, which are in physical connection with the object whose position beiiig 
monitored. Time based triangulation and (or) phase shift triangulation are used in such 
cases to periodically determine the position of the monitored object, radioactive 
emission probe 22 in this case. Examples of such a position ti^cking systems 
en^Ioyed in a variety of contexts using acoustic* (e.g., ultrasound) electromagnetic 

25 radiation (e.g., infrared^ radio frequency) or magnetic field and optical decoding are 
disclosed in, for example, U.S. Pat Nos. 5,412,619; 6,083,170; 6,063,022; 5,954,665; 
5,840,025; 5,718,241; 5,713,946; 5,694,945; 5,568,809; 5,546,951; 5,480,422 and 
5,391,199, which are incorporated by reference as if fiilly set forth herein. 

Radioactive emission probes are well known in the art and may use any one of 

30 a number of approaches for the determination of the amount of radioactive emission 
emanating from an object or portion thereof Dqiending on the type of radiation, such 
detectors typically include substances which when interacting with radioactive decay 
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emitted particles emit either electrons or photons in a level which is proportional over 
a wide linear range of operation to the level of radiation impinging thereon. The 
emission of electrons or photons is measurable and therefore serves to quantitatively 
determine radiation levels. Solid-state detectors in the form of N-type, P-type, 
5 PIN-type pixellated or unpixellated include^ for example, Ge, Si, Cdfe, CdZnTc, 
CdSe, CdZnSe, Higl2, TlBrl, GaAs, Inl, GaSe, Diamond, TlBr, Pbfc. hiP, ZnTe, HgBrl, 
a-Si, a-Sc, BP, GaP, CdS, SiC, AlSb, PbO, Bil3 and ZnSe detectors. Gas (e.g., CO2 
CH4) filled detectors include ionization chamber detectors, proportional chamber 
detectors and geiger chamber detectors. Scintillation detectors include organic 

10 scintillators crystals and liquids, such as CuHio j Ci4H]2 , C]oHs, etc.. Plastics, 
NE102A, NEI04, NEl 10, Pilot U and inorganic scintillalors, such as Nal, Csl, BGO, 
LSO, YSO, BaF, ZnS, ZnO, CaW04 and CdW04. Also known are scintillation fiber 
detectors. Scintillator coupling include pholomultiptier tube (PMT) of the following 
types: side-on type, head-on type, hemispherical type, position sensitive type, 

15 icrochannel plate-photomultiplier (MCP-PMTs) and electron multipliers, or 
photodiodes (and photodiodes arrays), such as Si photodiodes. Si PIN photodiodes. Si 
APD, GaAs(P) photodiodes, GaP and CCD. 

Figure S shows a narrow angle or wide angle radioactive emission probe 22\ 
Narrow or wide angle radioactive emission probe 22* includes a narrow slit 

20 (collimator) so as to allow only radiation arriving from a predetermined angular 
direction (e.g., 1 - 280 wide angle, preferably 1 ° - 80 - nanow angle) to enter the 
detector. Narrow or wide angle radioactive emission probes especially suitable for the 
configuration shown in Figure 10 are manufactured, for example, by Neoprobe, 
Dublin, Ohio (www.neoprobe.com), USA, Nuclear Fields, USA (www.nufi.com) 

25 IntraMedical Imaging, Los Angeles, CA, USA (www.gammaprobe.com). 

As is shown in Figure 6, such a detector is typically used to measure 
radioactivity, point by point, by scanning over the surface of a radioactive object from 
a plurality of directions and distances. In the example shown, scans from four 
different directions are employed. It will be ^leciated that if sufficient radioactivity 

30 records are collected from different angles and distances, and the orientation and 
position in space of detector 22' is simultaneously mionitored and recorded during such 
scans, a three-dimensional model of a radioactive region can be reconstimted and its 
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position in space detennined. If two or more detectors are co-employed, as shown in 
the configuration of Figure 10, the results maybe collected faster. 

Figure 7 shows another example of a radioactive emission probe, a spatially 
sensitive (pixelated) radioactive emission probe 22- (such as a gamma camera). 

5 Detector 22", in effect, includes an array of multitude narrow angle detector units 23. 
Such an airangement is used in accordance with the teachings of the present invention 
to reduce the amount of measurements and angles necessary to acquire sufficient data 
so as to reconstitute a three-dimensional model of the radioactive object. Examples of 
spatially sensitive radioactive emission probes employed in a variety of contexts are 

!0 disclosed in, for example, U.S. Pat. Nos. 4,019,057; 4,550,250; 4,831,262; and 
5,521,373; which are incorporated by reference as if set forth herein. An additional 
example is the COMPTON detector flittp://www.ucl.ac,uk/MedPhvs/posters/ 
giulia/giulia.htm) . Figure 8 shows a scan optionally made by q>atially sensitive 
radioactive emission probe 22" (such as a gamma camera). 

15 A radioactive emission detector of particular advantages for use in context of 

the present invention is the Compton gamma probe, since, in the Compton gamma 
probe, spatial resolution is independent of sensitivity and it appears possible to exceed 
the noise equivalent sensitivity of coUimated imaging systems especially for systems 
with high spatial resolution. The Compton probe is a novel type of gamma-probe that 

20 makes use of the icinematics of Compton scattering to construct a source image 
without the aid of mechanical collimators. Compton imaging telescopes were first 
built in the 1970s for astronomical observations [V. Schoenfelder et al., Astrophysical 
Journal 217 (1977) 306]. The first medical imaging laboratory instrument was 
proposed in the early 1980s [M. Singh, Med. Phys. 10 (1983) 421]. The potential 

25 advantages of the Compton gamma probe include higher efficiency, 3-D imaging 
without detector motion, and more compact and lightweight system. In the Compton 
gamma probe, high-energy gamma rays are scattered fiom a first detector layer (or 
detectors array) into a second detector layer array. For each gamma, the deposited 
energy is measured in both detectors. Using a line drawn between these two 

30 detectors, the Compton scattering equation can be solved to determine the cone of 
possible direction about this axis on which the gamma ray must have entered the first 
detector. The intersection of cones from many events is then developed to locate 
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gamma ray sources in the probe's field-of-view. Obviously only coincident events aie 
considered, and the more accurately their energy can be detennined, the less 
uncertainty fliere is in the spatial angle of the arrival cone. The probe's electronic 
system is combining coincidence measurements across many detectors and detectors 

5 layers with a very good energy resolution. The choice of the geometry and the 
material of the first layer detector plays a major role in the system imaging capability 
and depends on (i) material eCficiency of single Compton events, in relation to other 
interactions; (ii) detector energy resolution; and (iii) detector position resolution. In 
paiticular, the overall angular resolution results from the combination of two 

10 components, related to the energy resolution and to the pixel volume of the detector 

Thus, as now afforded by the present invention, connecting a radioactive 
emission probe to a position tracking system, permits simultaneous radioactivity 
detecting and position tracking at the same time. This enables the accurate calculation 
of the shape, size and contour of the radiating object and its precise position in a 

15 system-of-coordinates. 

The present invention thus provides a method for defining a position of a 
radioactivity emitting source in a system-of-coordinates. The method is effected by (a) 
providing a radioactive emission probe which is in communication with a position 
tracking system; and (b) monitoring radioactivity emitted firom the radioactivity 

20 emitting source, while at the same time, monitoring the position of radioactive 
emission probe in the system-of-coordinates, thereby defining the image of the 
radioactivity emitting source in the system-of-coordinates. 

It will be s^preciated by one of skills in the art that the model produced by 
system 20 is projectable onto any of the other systems-of-coordinates, or alternatively, 

25 the system-of-coordinates defined by position tracking system 24 may be shared by 
other position tracking systems, as is further detailed hereinbelow, such that no such 
projection is required. 

Thus, as is fiirthcr shown in Figure 1, system 20 of the present inv^tion can be 
used for calculating a position of a radioactivity emitting source in a first 

30 system-of-coordinates 28 and fiirther for projecting the image of the radioactivity 
emitting source onto a second system-of-coordinates 28*. The system includes 
radioactive emission probe 22, position tracking system 24 which is connected to and 
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(or) communicating with radioactive emission probe 22, and data processor 26 which 
is designed and configured for (i) receiving data inputs from position tracking system 
24 and J&om radioactive emission probe 22; (ii) constructing flie image of the 
radioactivity emitting source in the first system-of-coordinates; and (iii) projecting the 

5 image of the radioactivity emitting source onto the second system-of-coordinates. 

A method for calculating a position of a radioactivity emitting source in a first 
system-of-coordinates and for projecting the image of ^e radioactivity emitting source 
onto a second system-of-coordinates is also offered by the present invention. This 
method is effected by (a) providing a radioactive emission probe being in 

10 communication with a position tracking system; and (b) monitoring radioactivity being 
emitted fiom the radioactivity emitting source, while at the same time, monitoring the 
position of the radioactive mission probe in the first system-of-coordinates, thereby 
defining the image of the radioactivity emitting source in the first 
system-of-coordinates and projecting the image of the radioactivity emitting source 

15 onto the second system-of-coordinates. 

It will be appreciated that the combination of a radioactive emission probe and 
a position tracking system connected thereto and (or) communicating therewith allows 
a suitable data processor to generate a two- or three-dimensional image of the 
radioactivity emitting source. An algorithm can be used to calculate image intensity 

20 based on, for example^ a probability function which averages radiation counts and 
generates an image in which the shorter the time interval between radioactive counts, 
the brighter the image and vise versa, while down-compensating when a location is 
re-scanned» A fi*ee-hand scanning with a directional detector can be employed for this 
purpose. 

25 hi one embodiment, when scanning a body area with the detector, the detector 

is made to follow a three-dimensional surface, which defmes the body curvature and in 
effect is used also as a position tracking pointer. This information can be used to 
define the position of the radioactive source with respect to the outer surface of the 
body, so as to create a three-dimensional map of both the radioactive source and of the 

30 body curvature. This approach can also be undertaken in open surgeries, such as open 
chest surgeries so as to provide the surgeon in real time with information concerning 
the functionality of a tissue. 
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The radioactive emission probe, which can be used in context of the present 

invention can be a beta emission detector, a ganrnia omission detector, a positron 

emission detector or any combination thereof. A detector that is sensitive to both beta 

and (or) positron and gamma emission can be used to improve localization by sensing 

S for example gamma mission distant from the source and sensing beta or positrons 

emission closer to the source. A beta detector is dedicated for the detection of either 

131 

electrons from sources such as Iodine^ or positrons from sources such as 
18 

Fluorine. A gamma detector can be designed as a single ^ergy detector or as a 
detector that can distinguish between different types of energies, using the light 

10 intensity in the scintilator as a relative measure of the gamma energy. Also, the 
detector can be designed to utihze coincidence detection by using detectors facing one 
another (180 degrees) with the examined organ or tissue in-between. The radiation 
detector can have different collimators with different diameters. A large bore will be 
used for high sensitivity with lower resolution while a small bore collimator wH have 

15 higher resolution at the expense of lower sensitivity. 

Another possibility is to have a the collimator moving or rotating with the 
opening eccentric so that a different solid angle is exposed to the incoming photons at 
any one time, thus gathering the photons from overlapping volumes at different time 
intervals. The rest of the image processing is similar if the probe moves or if the 

20 collimator eccentric opening moves. 

System 20 of the present invention can be used in concert with other medical 
devices, such as, but not limited to, any one of a variety of imaging modalities and (or) 
surgical instruments. 

Structural imaging modalities, which provide anatomic, or structural maps of 

25 the body, are well known in the art. The main modalities that serve for two- 
(projectional or cross sectional) or three- (consecutive cross sectional) dimensional 
imaging are a planer X-ray imager, a fluoroscope, a computerized tomography scanner, 
a magnetic resonance imager an ultrasound imager, an impedance imager, and an 
optical camera. 

30 Medical images taken of the human body are typically acquired or displayed in 

three main orientations (i) coronal orientation: in a cross section (plane), for example, 
across the shoulders, dividing the body into front and back halves; (ii) sagittal 
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orientation: in a cross section (plane), for example, down the middle, dividing the body 
into left and right halves; and (iii) axial orientation: in a cross section (plane), 
perpendicular to the long axis of the body, dividing the body into upp^ and lower 
halves. Oblique views can also be acquired and displayed. 
5 Various types of X-ray imaging are central to diagnosis of many types of 

cancer. Conventional X-ray imaging has evolved over the past 1 00 years, but the basic 
principal is still the same as in 1895, when first introduced. An X-ray source is turned 
on and X-rays are radiated through the body part of interest and onto a film cassette 
positioned under or behind the body part. The energy and wavelength of the X-rays 

10 allows them to pass througih the body part and create the image of the internal 
structures like bones. As the X-rays pass through the hand, for instance, they are 
attenuated by the different density tissues they encounter. Bone attenuates a great deal 
more of the X-rays than the soft tissue surrounding it because of its grater density. It is 
these differences in absorption and the corresponding varying exposure level of the 

15 film that creates the images. In fact. X-ray imaging results in a projection of the 
integrated density of column-voxels defined by the X-rays as they pass through the 
body. 

Fluoroscopy is a method based on the principals of film X-ray that is useful for 
detecting disorders and tumors in the upper gastro-inlestinal (GI) system (for example, 

20 the stomach and intestines). Fluoroscopic imaging yields a moving X-ray picture. The 
physician can watch the screen and see an image of the patient's body (for example the 
beating heart). Fluoroscopic technology improved greatly with the addition of 
television cameras and fluoroscopic "image intensifiers**. Today, many conventional 
X-ray systems have the ability to switch back and forth between the radiographic and 

25 fluoroscopic modes. The latest X-ray systems have the ability to acquire the 
radiograph or fluoroscopic movie using digital acquisition. 

Computed Tomography (CT) is based on the X-ray principal, where the fikn is 
replaced by a detector that measures the X-ray profile. Inside the covers of the CT 
scanner is a rotating frame which has an X-ray tube mounted on one side and the 

30 detector mounted on the opposite side. A tan beam of X-ray is created as the rotating 
frame spins the X-ray tube and detector around the patient. Each time the X-ray tube 
and detector make a 360° rotation, an image or "slice" has been acquired. This "slice" 
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is coUimated to a thickness between 1 mm and 10 mm using lead shutters in front of 
the X-ray tube and X-ray detector. 

As the X-ray tube and detector make this 360"^ rotation, the detector takes 
num^us i^files of the attenuated X-ray beam. Typically, in one 360^ lap, about 

5 1,000 profiles are sampled. Each profile is subdivided spatially by the detectors and 
fed into about 700 individual channels. Each profile is then backwards reconstructed 
(or "back projected") by a dedicated computer into a two-dimensional image of the 
"slice" that was scanned. 

The CT gantiy and table have multiple microprocessors that control the 

10 rotation of the gantry, movement of the table (up/down and in/out), tilting of the gantry 
for angled images, and other functions such as turning the X-ray beam on an off. The 
CT contains a slip ring that allows electric power to be transferred from a stationary 
power source onto the continuously rotating gantry. The innovation of the power slip 
ring has created a renaissance in CT called spiral or helical scanning. These spiral CT 

15 scanners can now image entire anatomic regions like the lungs in a quick 20 to 30 
second breath hold. Instead of acquiring a stack of individual shces which may be 
misaUgned due to slight patient motion or breathing (and lung/abdomen motion) in 
between each slice acquisition, spiral CT acquires a volume of data with the patient 
anatomy all in one position. This volume data set can then be computer-reconstructed 

20 to provide three-dimensional models such as of complex blood vessels like the renal 
arteries or aorta. Spiral CT allows the acquisition of CT data that is perfectly suited 
for three-dimensional reconstruction. 

MR hnaging is superior to CT in detecting soft tissue lesions such as tumors as 
it has excellent contrast resolution, meaning it can show subtle soft-tissue changes 

25 with exceptional clarity. Thus, MR is often the method of choice for diagnosing 
tumors and for searching for metastases. MR uses magnetic energy and radio waves to 
create single or consecutive cross-sectional images or *'slices" of the human body. The 
main component of most MR systems is a large tube shaped or cylindrical magnet 
Also, there are MR systems with a C-sliaped magnet or other type of open designs. 

30 The strength of the MR systems magnetic field is measured in metric units called 
"Tesla". Most of the. cylindrical magnets have a strength between 0.5 and 1.5 Tesla 
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and most of the open or C-shaped magnets have a magnetic strength between 0.01 and 
0.35 Tesla. 

Inside the MR system a magnetic field is created. Each total MR examination 
typically is comprised of a series of 2 to 6 sequences. An "MR sequence" is an 

5 acquisition of data that yields a specific image orientation and a specific type of image 
appearance or "contrast". During the examination, a radio signa] is turned on and off, 
and subsequently the energy which is absorbed by different atoms in the body is 
echoed or reflected back out of the body. These echoes are continuously measured by 
"gradient coils" tliat are switched on and off to measure the MR signal reflecting back. 

10 In the rotating frame of reference, the net magnetization vector rotate from a 
longitudinal position a distance proportional to the time length of the radio frequency 
pulse. After a certain length of time, the net magnetization vector rotates 90 degrees 
and lies in the transverse or x-y plane. It is in this position that the net magnetization 
can be detected on MRI. The angle that the net magnetization vector rotates is 

13 commonly called the 'flip' or tip* angle. At angles greater than or less than 90 degrees 
tfaoie will still be a small component of the magnetization that will be in the x-y plane, 
and therefore be detected. Radio frequency coils are the "antenna" of the MRI system 
that broadcasts the RF signal to the patient and (or) receives the return signal. RF coils 
can be receive-only, in which case the body coil is used as a transmitter; or tran^t 

20 and receive (transceiver). Surface coils are the simplest design of coil. They are 
simply a loop of wire, either circular or rectangular, that is placed over the region of 
int^est. 

A digital computer reconstructs these echoes into images of the body. A 
benefit of MRI is that it can easily acquire direct views of the body in almost any 

25 orientation, while CT scanners typically acquire cross-sectional images perpendicular 
or nearly perpendicular to the long body axis. 

Ultrasound imaging is a versatile scanning technique that uses sound waves to 
create images of organs or anatomical structures in order to make a diagnosis. The 
ultrasound process involves placing a . small device called a transducer, against the skin 

30 of the patient near the region of interest, for example, against the back to image the 
kidneys. The ultrasound transducer combines functions of emitting and receiving 
sound. This transducer produces a stream of inaudible, high frequency sound waves 
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which p^etrate into the body and echo off the organs inside. The transducer detects 
sound waves as they echo back fiom the internal structures and contours of the organs. 
Different tissues reflect these sound waves differently, causing a signature which can 
be measured and transformed into an image. These waves are received by the 
S ultrasound machine and turned into live pictures with the use of computers and 
reconstruction software. 

Ultrasound scanning has many uses, including: diagnosis of disease and 
structural abnormalities, helping to conduct other diagnostic procedures, such as 
needle biopsies etc. 

10 There are limitations to some ultrasound techniques: good images may not be 

obtained in every case, and the scan may not produce as precise results as some other 
diagnostic imaging procedures. In addition, scan results may be affected by physical 
abnormalities^ chronic disease, excessive movement, or incorrect transducer 
placement. 

15 Both two- (cross sectional) and three- (consecutive cross-sectional) ultrasound 

imaging techniques are available nowadays. Worth mentioning is the Doppler 

three-dimensional ultrasound imaging. 

In many cases imaging modalities either inherently include (e.g., fluoroscopy 

CT, MRT) and (or) are integrated with position-tracking-systems, which enable the use 
20 of such systems to reconstruct three-dimensional image models and provide their 

position in a system-of-coordinates. 

It will be appreciated that, similar to the vision system, also an optical camera 

can be used to generate three-dimensional imagery date according to the present 

invention by imaging a body from a plurality (at least two) directions. This type of 
25 imaging is especially applicable in open chest surgeries or other open surgeries. 

Software for calculating a three-dimensional image from a pair of stereoscopic images 

is well known in the art. 

Thus, as used herein and in the claims section that follows, the phrase 

''three-dimensional imaging modality'* refers to any type of imaging equipment which 
30 includes software and hardware for generating a three-dimensional image. Such an 

equipment can generate a three-dimensional image by imaging successive 

cross-sections of a body, e.g., as if viewed from a single direction. Alternatively, such 
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an equipment can generate a three-dimensional image by imaging successive 
cross-sections of a body, e.g., as if viewed ftom a single direction. Alternatively, such an 
equipment can generate a three-dimensional image by imaging a body from different 
angles or directions (typically two angles) and thereafter combining the data into a 

5 three-dimensiona) image. 

In accordance with the present invention, a structural hnaging probe, for example, 
an ultrasound probe, may be combined with a position tracking system, in a manner 
analogous to System 20 of radioactivity emission detector 22 and position tracking system 
24 (Figure I), Smiilarly, a miniature MRI probe, for example, as taught by U.S. Patent 

10 5,572,132, to Pulyer, et al., entitled, ''MRI probe for external imaging," whose disclosure 
is incorporated herein by reference, which teaches an MRI catheter for endoscopic 
imaging of tissue of the artery wall, rectum, urinal tract, intestine, esophagus, nasal 
passages, vagina and other biomedical applications, may be used. 

Additionally, a structural imaging probe, for example, an ultrasound probe, may 

15 be combined with system 20 of radioactivity emission detector 22 and position tracking 
system 24 (Figure 1), so as to have: 

• structural imaging; 

• fimctional imaging; and 

• position tracking, 

20 simultaneously. Similarly, other structural imaging modalities, as known, may be 
combined with system 20 of Figure 1. 

Furthermore, the structural imaging modality may have an additional, dedicated 
position tracking system, or share position tracking system 24 of radioactivity emission 
detector 22. 

25 Position tracking may also be accomplished by using imaging information from 

the structural imaging system. This may be accomplished by tracking relative changes 
from one image to another to determine relative motion. 

Surgical instruments are also well known in the art and may use any one of a 
plurality of configurations in order to perform minimally-invasive surgical procedures. 

30 Examples include laser probes, cardiac and angioplastic catheters, endoscopic probes, 
biopsy needles, aspiration tubes or needles, resecting devices, ultrasonic probes, fiber 
optic scopes, laparoscopic prob s, thermal probes and suction/irrigation probes. 
Examples of such surgical instruments employed in a variety of medical cont xts are 
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disclosed in, for example, U.S. Pat. Nos. 6,083,170; 6,063,022; 5,954,665; 5,840,025; 
5,718,241; 5,713,946; 5,694,945; 5.568,809; 5,546,951; 5,480,422 5,391,199, 
5,800,414; 5,843,017; 6,086,554; 5,766,234; 5,868,739; 5,911,719; 5,993,408; 
6.007,497; 6,021,341; 6.066,151; 6,071,281; 6,083,166 and 5,746,738, which are 
3 incorporated by reference as if fiilly set forth herein. 

For some applications^ examples of which are provided in the list of patents 
above, surgical instruments are integrated with position-tracking-systems, which 
enable to monitor the position of such instruments while placed in and guided through 
the body of a treated patient. 

10 According to a preferred embodiment of the present invention, the surgical 

instrument is equipped with an additional radioactive emission probe attached thereto 
or placed therein. This additional detector is used, according to preferred 
embodiments of the invention^ to fine tune the location of radioactive emission fix)m 
Avithin the body, and in closer proximity to the radioactive source. Since the surreal 

1 5 tool is preferably in communication with a position-tracking system, the position of the 
additional detector can be monitored and its readouts used to fine tune the position of 
the radioactive source within the body. Thus, according to this aspect of the present 
invention, at least one extracorporeal detector and an intracorporeabeal detector are 
used in concert to determine the position of a radioactive source in the body in hig^st 

20 precision. The extracorporeal detector provides the general position of the source and 
is used for directing the surgical instrxmient thereto, whereas the Intracorporeal 
detector is used for reassuring prior to application of treatment or retrieval of biopsy 
that indeed the source was correctly targeted at the highest precision. 

While according to a presently preferred embodiment of the invention two 

25 detectors, one extracorporeal md one Intracorporeal, are employed as described above, 
for some applications a single Intracorporeal detector may be employed, which 
detector is attached to or integrated with a surgjcal instrum^t whose position is 
tracked. 

The use of Intracorporeal and extracorporeal detectors calls for careful choice 
30 of the radioactive isotope employed with the radiopharmaceutical. While the 
extracorporeal detector can be constructed with a suitable collimator for handling 
strong radiation, such as ganuna radiation, the Intracorporeal detector is miniature by 
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nature and is limited in design and construction by the construction of the surgical 

instrument with which it is employed. Since collimators for high energy (80 - 5 1 1 

KeV) gamma radiation are robust in nature, they are not readily eogageable with 

miniature detectors. Electron (beta) and positron radiation are characterized by: (i) 

5 they highly absorbed by biological tissue as they are of lower energy and higher 

chemical reactivity; and (ii) they are readily collimated and focused by thin metal 

colhmators. It is also possible to use low energy gamma radiation (10-30 KeV) for 

Intracorporeal applications since the collimation of these gamma photons can be 

achieved with thin layers of Tantalum or Tungsten. As such, the radio phannaceutical 

10 of choice is selected to emit both gamma and beta and (or) positron radiation, whereas 

the extnicoiporeal detector is set to detect the high energy gamma radiation, whereas 

the Intracorporeal detector is set to detect the low energy gamma, beta and (or) 

positron radiation. Isotopes thai emit both high eneigy gamma and (or) low energy 

gamma, beta and (or) positron radiation and which can be used per se or as a part of a 

18 111 123 

15 compound as radiopharmaceuticals include, without limitation, F, In and I in 

18 

radiopharmaceuticals, such as, but not limited to, 2-[ F]fluon>-2-deoxy-D-glucose 
(^^FDG), ^ * ^In-Pentetreotide ([^ ^ ^In-DTPA-D-Phe^]^treotide), 

L-3-[ I]-Iodo-alpha-methyl-tyrosine (IMT), 0-(2-[ F]fluoroethyl)-L-tyrosine 
(L-(^^F]FET), ^^^hi'Capromab Pendetide (CYT-356. Prostascint) and 

20 ^ ^ ^ In-Satumomab Pendetide (Oncoscint). 

Figure 1 1 illustrates a system in accordance with this aspect of the present 
invention. A surgical instrument 100 is shown connected to a resection/aspiration 
control element 102 as well known in the art. Surgical instrument 100 includes a 
radioactive emission probe 104, which has a collimator 106 for collimating low energy 

25 gamma, beta and (or) positron radiation. In some embodiments, as indicated by airow 
108, detector 104 may be translated within instrument 100. A position tracking system 
having one elematit thereof 110 attached to instrument 100 and another element 
thereof 1 12 at a iixed location serves to monitor the position of instrument 100 at all 
times in two, three and up to six degrees of freedom. Radioactive emission probe 104 

30 communic^ates with a counter 114 for counting low energy gamma, beta and (or) 
positron radiation. All the data is communicated to, and processed by, a processor 
1 1 6. The 2D or 3D data may be projected and displayed along with 2D or 3D imaging 
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data derived from an imaging modality using a shared presentation device as described 
elsewhere herein. A real or virtual image of the surgical instrument itself may also be 
co-displayed. Examples of commercially available radiation emission detectors that 
can fit inside, for example, a biopsy needle include scintillating plastic optical fibers 
5 like SlOl and S104, manufactured by PPLASTIFO or an optical fiber communicating 
with a scintillator (either detector paint or scintillation crystal) at the fiber edge. The 
level of detected radiation can be reported visually or by an audio signal, as is well 
known in the art. 

Thus, a surgical instrument equipped with a radiation emission detector and 
10 which is connected to and (or) communicating with a position tracking system forms 
one embodiment of this aspect of the present invention. Such a design acting in 
concert with either conventional imaging modalities and (or) extracorporeal radiation 
emission detectors form other embodiments of this aspect of the invention* In all 
cases, a surgical instrument equipped with a radiation emission detector and which is 
15 connected to and (or) communicating with a position tracking system serves for in situ 
fine tuning of a radioactive source in the body. 

It will be appreciated that in some minimally-invasive procedures even the 
position of the patient him or herself is monitored via a position tracking system, 
using, for example, electronic or physical fiducial markers attached at certain locations 
20 to the patient* s body. 

Thus, as is fiirther detailed hereinbelow, by projecting the three-dimensional 
data and positions received bom my of the above mentioned devices into a common 
system of coordinates, or alternatively, employing a common position tracking system 
for all of these devices, one can integrate the data into a far superior and 
25 comprehensive presentation. 

An example to this desired outcome is shown in Figure 9. In the embodiment 
shown, four indq>endent position tracking systems 50, 52, 54 and 56 are used to track 
the positions of a patient 58, an imaging modality 60, a radioactive emission probe 62 
and a surgical instrument 64 in four independent sy$tems-of-coordinates 66, 68, 70 and 
30 72, respectively. If the patient is still, no tracking of the patient's position is required. 

It will be appreciated that any subset or all of the position tracking systems 
employed may be integrated into one or more common position tracking systems, and 
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(or) that any subset or all of the position tracking systems employed may share one or 
more systemE-of-coordinates» and fiirther that any positional data obtained by any of 
the position tracking systems described in any of the systems-of coordinates may be 
projected to any other system of coordinates or to an independent (fifth) system of 

5 coordinates 74. Tn one preferred embodiment, applicable for applications at the torso 
of the patient, the system of coordinates is a dynamic system of coordinates which 
takes into account the chest breathing movements of the patient during the procedure. 

As indicated at 76, the raw data collected by detector 62 is recorded and, as 
indicated at 78, the position and the radioactive data records are used to generate a 

10 three-dimensionai model of a radiophaimaceuticai uptake portion of a body 
component of the patient. 

Similarly, as indicated at 80, the imagery data collected by imaging modality 
60 is recorded and the position and the imagery data records are used to generate a 
three-dimensional model of the imaged body component of the patient. 

!5 All the data collected is then fed into a data processor 82 which processes the 

data and, as indicated at 84, generates a combined or superimposed presentation of the 
radioactive data and the imagery data, which is in positional context with patient 58 
and surgical instrument 64. 

Instrument 64, which by itself can be presented in context of the combined 

20 presentation, may then be used to perform the procedure most accurately. Processor 
82 may be a single entity or may include a plurality of data processing stations which 
directly communicate with, or even integral to, any one or more of the devices 
described. 

Additionally or alternatively, a structural imaging probe, for example, an 
25 ultrasound probe, or another structural probe, as known, may be incorporated with the 
surgical instrument. 

The present invention provides a major advantage over prior art designs 
because it positionally integrates data pertaining to a body portion as retrieved by two 
independent imaging techniques, conventional imaging and radioactive imaging, to 
30 thereby provide a surgeon with the ability the fine point the portion of the body to be 
sampled or treated. 
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It will be appreciated that subsets of the devices described in Figure 9 may be 
used as stand-alone systems. For example, a combination of detector 62 with its 
position-tracking system and instrument 64 with its position-tracking-system may in 
some instances be sufficient to p^form intrabody procedures. For mere diagnostic 
5 purposes^ without biopsy^ a combination of detector 62 position-tracking-system and 
* modality 60 position-tracking-system are sufficient 

Reference is now made to Figure 12, which illustrates an imaging system 200 
constracted and operative in accordance with a preferred embodiment of the present 
invention, baagtng system 200 preferably includes a radiation probe 202, such as the 
10 narrow angle radioactive emission probe 22' described hereinabove with reference to 
Figures 5 and 10. 

A position sensor 204 is provided for sensing the position of radiation probe 
202. Position sensor 204 may be physically attached to radiation probe 202, or may be 
distanced therefipom. Position sensor 204 transmits the sensed position data to a 
15 position tracking system 206. Position tracking system 206 may be a system like 
position tracking system 24, described hereinabove with reference to Figure 1, and 
position sensor 204 may be any kind of sensor applicable for such position tracking 
systems. 

Another method which can be used to locate the source of radiation emission is 
20 by using a small hand held gamma camera 205 (such as the DigiRad 2020tc Imager 
TM, 9350 Trade Place, San Diego, California 92126-6334, USA), attached to position 
sensor 204. 

Position tracking system 206 enables radiation probe 202 to freely scan back 
and forth in two- or three-dimensions over the area of interest of the patient, preferably 
25 incrementing a short distance between each scan pass. Position tracking system 206 
tracks flie position of radiation probe 202 with respect to a position tracking coordinate 
system, such as Xp, Yp and Zp^ with an origui Op. 

Imaging system 200 also includes a medical imaging system 208, such as, but 
not limited to, computed or computerized tomography (CT), magnetic resonance 
30 imaging (MRl), ultrasound imaging, positron emission tomography (PET) and single 
photon emission computed tomography (SPECT), for example. Medical imaging 
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system 208 provides images of a patient 209 with respect to a medical imaging 
coordinate system, such as Xm, ^ origin Om. 

Imaging system 200 also includes a coordinate registration system 210, such as 
that described in commonly owned, now abandoned, U.S. Patent Application No. 
S 09/610^490, entitled, ''method for registering coordinate systems," the disclosure of 
which is incoiporated herein by reference. Coordinate registration system 210 is 
adapted to register the coordinates of the position traddng coordinate system with 
those of the medical imaging coordinate system. 

Position tracking system 206, medical imaging system 208 and coordinate 

10 registration system 210 are preferably in wired or wireless connnunication with a 
processing unit 2 1 2 (also referred to as a data processor 2 1 2). 

In operation of imaging system 200, after administration of a 
radiopharmaceutical to patient 209, a clinician/physician/surgeon (not shown) may 
move or scan radiation probe 202 about a target area under examination. A 

15 physiological activity map of the target area is obtained by measuring the radiation 
count rate with radiation probe 202, and by correlating the count rate with the count 
rate direction with position tracking system 206, which follows the motion of the 
moving or scanning radiation probe 202. 

Refo^ence is now made to Figure 13, which illustrates image formation with 

20 radiation probe 202, in accordance with a preferred embodiment of the present 
invention. For the purposes of simplicity, the example shown in Figure 13 is for a 
single dimension image formation, but it is readily understood that the same prmciples 
hold true for any other dimensional image formation. 

In one exanq^le of carrying out the invention, radiation probe 202 may be a 

25 gamma ray detector probe that comprises a collimator 211 and radiation detector 213. 
Gamma rays are projected through the probe collimator 211 onto radiation detector 
213, which produces electronic signals in accordance vtdth the radiation detected. 
Radiation probe 202 sends pulses to a probe counter 215 which may include a pulse 
height analyzer circuit (not shown). The pulse height analyzer circuit analyzes the 

30 electronic signals produced by radiation detector 213. If the electronic signals are 
within a selected energy window, the level of radiation, i.e., number of radiation 
counts, is counted by probe counter 215. 



10/07 03 THU 13:33 FAX 972 3 6127575 



G. E. EHRLICH 



^036 



34 

Examples of suitable radiation detectors include a solid state detector (SSD) 
(CdZnTe, CdTe, Hgl, Si, Ge, and the like), a scintillation detector (NaI(TlX LSO, 
GSO, Csl, CaF, and the h'ke), a gas detector, or a scintillating fiber detector (SlOl^ 
SI 04, and the like), for example. 

5 The position sensor 204 associated with the radiation probe 202 senses the 

position of radiation probe 202, and position tracking system 206 calculates and 
monitors the motion of radiation probe 202 with respect to the position tracking 
coordinate system. The motion is calculated and monitored in two, three and up to six 
dimensions * the linear directions of the X, Y and Z axes as well as rotations about the 

10 X, Y and Z axes, i.e., rotational angles p, 9 and <p, req)ectively. 

Examples of suitable position tracking systems include a measurment 
mechanical aim (FaroAnn, http:/ywww.faro.com/products/faioann.asp), optical 
(racking systems (Northern Digital Inc., Ontario, Canada NDI-POLARIS passive or 
active systems), magnetic tracking systms (NDl-AURORA), injfrared tracking 

15 systems (E-PEN system, http://www.e-penxomX and ultrasonic tracking systems 
(E-PEN system), for example. 

Processing unit 212 combines the radiation probe count rate from probe 
counter 215 together with the positional infonnation from position tracking system 
206, and uses an imaging software algorithm 217 to form a two-dimoisional or 

20 three-dimensional radiotracer-spread image of the target area inside the patient's body. 
The ^atial probe positions together with the spatial count rates may be stored in 
memory or displayed on a computer monitor 214 as a pattern of marks corresponding 
to the spatial and count rate position. 

An example of such a pattern is shown in Figure 14, which illustrates a 

2S single-dimensional, unprocessed simulation of a radiation point source 218 (Figure 
13), 30 mm deep inside the human body, detected by using a 10 mm nuclear radiation 
probe 202 coupled to position tracking system 206. The graph of Figure 14 indicates 
to a physician that there is a peak count rate of about 500 in the probe position of about 
SO mm. 

30 In one embodiment of the inv^tion, the imaging software algorithm 217 

employs an averaging process to refme the curve of Figure 14. This averaging process 
will now be described with reference to Figure 15. 
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Probe counter 21 S feeds probe count rate information H(Xc, Yc, Zc, 9, q> ) to 

« 

processing xmit 212 (step 301). Position sensor 204 feeds probe position information 
(Xc, Yc, Zc, p, 0, <p) to processing unit 212 (step 302). Probe parameters (such as its 
physical size, dx, dy, dz) are also input into processing unit 212 (step 303). 

5 Processing unit 212 then finds all the voxels (i.e., volume pixels) that represent 

the probe volume in the processing unit memory (step 304), i.e., Xc+dx, Yc+dy, 
Zc+dz. Processing unit 212 calculates the number of times that the calculation process 
has been done in each voxel from the beginning of the image formation (step 305), i.e., 
M(Xc+dx, Yc+dy, Zc+d;^). Processing unit 212 then calculates the new average count 

10 rate values in each voxel (step 306), in accordance with the formula: 

N(Xc+dx, Yc+dy, Zc+dz) = 

[ N(Xc+dx, Yc+dy, Zc+dz) + N(Xc, Yc, Zc, p, 0, (p) ] / 

[ M(Xc+dx, Yc+dy, Zc4dz) +1 ] 
Processing unit 212 then corrects the display image that represents the 
15 perceived voxels at N(Xc+dx, Yc+dy, Zc+dz) (step 307). The algorithm then r^ats 
itself for the next probe position (step 308). 

llie resulting graph of the averaging algorithm of Figure 15, as applied to the 
example of Figure 14, is shown in Figure 16. 

Figures 17 and 18 respectively show examples of a hot cross phantom image 
20 and a hot 4.77 mm bar phantom image, produced by a gamma radiation probe coupled 
with position tracldng system 206 and the averaging algorithm of Figure 15. The 
probe images were made by using a probe, EG&G Ortec NaI(Tl) model 905-1 
(thickness^l", diameter =1") connected to a ScintiPack model 296. The position 
tracking system used was the Ascension miniBIRD, commercially available from 
25 Ascension Technology Corporation, P.O. Box 527, Burlington, Veraiont 05402 USA 
(http;//www.ascension-tech.com/graphic.htm). The magnetic tracking and location 
systems of Ascension Technology Corporation use DC magnetic fields to overcome 
blocking and distortion from nearby conductive metals. Signals pass through the 
human body without attenuation. 
30 In another embodiment of the invention, the imaging softwarc algorithm 217 

may employ a minimizing process to refine the curve of Figure 14 as is now described 
with reference to Figure 19. 
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Probe counter 215 feeds probe count rate infonnation ]M(Xc, Yc, Zc, p, 0, cp) 
processing unit 212 (step 401). Position sensor 204 feeds probe position infonnation 
(Xc, Yc, Zc, p. 9. <p) to processing unit 212 (step 402). Probe parameters (such as its 
physical size, dx, dy, dz) are also input into processing unit 212 (step 403). 

5 Processing unit 212 then finds all the voxels that represent the probe volume in 

the processing unit memory (step 404), i.e., Xc+dx, Yc+dy, Zc4dz. From the voxels 
that represent the probe volume in the processing unit memory, processing unit 212 
finds those that have a higher count rale value N(Xc+dx, Yc+dy, Zc+dz) than the 
inputted probe count rate N(Xc, Yc. Zc, p, 9, (p) (step 405). Processing unit 212 then 

10 changes the higher count rate voxels to that of inputted probe count rate N(Xc, Yc, Zc, 
p, G, <i>) (step 406), and corrects the display image at the higher count rate voxels 
N(Xc+dx, Yc+dy, Zc+dz) (step 407). The algorithm then repeats itself for the next 
probe position (step 408). 

The resulting graph of the minimizing algorithm of Figure 19, as applied to the 

15 example of Figure 14, is shown in Figure 20. 

Reference is now made to Figure 21, which illustrates an image reconstruction 
system 450, constructed and operative in accordance with a preferred embodiment of 
the present invention. Image reconstruction system 450 produces a combined image 
451 made up of the images coming fiom the medical imaging system 208 with the 

20 position of the peak radiation location (and its uncertainty area) fi^om processing imit 
212, together with the location of a therapeutic instrument 452, such as a biopsy 
needle. The combined image 451 allows the. physician to better assess the relative 
position of therapeutic instrument 452 in relation to the anatomical image (from 
medical imaging system 208) and the position of the radioactive area as inferred by the 

25 radiation detection algorithm. 

In accordance widi preferred embodiments of the present invention, image 
acquisition and reconstruction algorithms are provided, based on image acquisition 
with wide-aperture collimation and image reconstruction, which includes 
deconvolution. The overall algorithms are heiethereto referred to as wide-aperture 

30 collimation - deconvolution algorithms. 

In essence, the wide-;^ture collimation - dec<Hivolution algorithms enable one 
to obtain a hi^i-efficiency, high resolution image of a radioactivity emitting source, by 



10/07 03 THU 13:35 FAX 972 3 6127575 



G. £. £HRLICH 



@039 



37 

scanning the radioactivity emitting source with a probe of a wide-aperture collimator, 
and at the same time, monitoring the position of the radioactive emission probe, at 
vay fine time intervals, to obtain the equivalence of (ine-aperture collimation. The 
bluiring effect of flie wide aperture is then corrected mathematically. 

5 The wide-apeture collimation • deconvolution algorithms are described 

hereinbelow, in conjunction with Figures 27A - 271, and Figure 22. Experimental 
results, showing images produced by the mdc-Hpctare collimation - deconvolution 
algorithms of the present invention, in comparison to a conventional gamma camera 
are seen in Figures 23 A - 26B, and Figure 36, hereinbelow. 

10 It will be appreciated that the wide-apeture collimation - deconvohition 

algorithms of the present invention apply also to systems on fixed gantries, such as 
gamma camera and SPECT. 

Referring further to the drawings. Figures 27A-27I describe wide-aperture 
collimation - deconvolution algorithms, in accordance witli preferred embodiments of 

15 the present invention. The wide-ap^ture collimation - deconvolution algorithms are 
designed for estimating the distribution of radiation sources in a control volume, thus 
constructing an image of the radiation sources in the control volume. For simplicity, it 
is assumed that the radiation sources comprise dot sources that radiate uniformly in all 
directions, are localized, and arc smoothly distributed in the control volume, 

20 Consider a single-pixel, wide-bore-collimator probe 707, seen in Figure 27H, 

formed for example, as a lead tube collimator 708, having a length, LcoUhmior of about 
20 mm and a radiation detector 706, for example, a solid state CdZnTI crystal, at its 
base, with a diameter, D^eM:\xx of about 10 mm, so that the ratio of Lcoiimiator to Ddeuxiw 
is about 2, An angular opening 8 of wide-bore-collimator probe 707 may be about 40 

25 degrees; wide-bore-collimator probe 707 is sensitive to any incident photons within 
the confining of about a 40-degree angular opening. It will be appreciated that 
single-pixel, wide-bore-colUmator probes of other dimensions and ratios of Lcoiiiimtor 
to Ddttcctor are similarly possible. For example, the ratio of Lconm»ter to Ddctector niay be 
about 1, or about 3. 

30 By contrast, as seen in Figure 271, conventional gamma cameras have a 

multi-pixel, small-aperture collimation probe, wherein a collimator - pixel 
arrangement 709 is ahout 2 nun in diameter and 40 mm in length, so that the ratio of 
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Lcojiinaior to Dtfctector is about 20. An angular opening e of the small-aperture 
collimation probe is very small. In consequence, wide-bore-collimator probe 707 may 
be about 30 times more sensitive than the small-aperture coUimation probe of 
collimator - pixel arrangement 709. 

« 

5 Yet, there are disadvantages to wide-bore-collimator probe 707: 

1 . a single-pixel probe does not lend itself to the generation of an image; and 

2. a wide aperture bluis the information regarding the direction from which the 
radiation comes. 

With regard to the first point, the generation of an image, as 
10 wide-bore-collimator probe 707 is moved in space, while its position and angular 
orientation are accurately tracked, at very short time intervals, for example, of about 
100 or 200 ms, or any other short time interval, a one-, two-, or three-dimensional 
image of counting rate as a ftmction of position can be obtained, by data processing. 
More specifically, knowing the position and orientation of probe 707 at each 
15 time interval and the photon count rate at that position and orientation can be used to 
reconstruct the radiation density of an unknown object, even in three dimensions. 

With regard to the second point, the blurring effect of the wide-angle aperture, 
while it is known that the infonnation is blurred, or convolved, by the wide-aperture 
collimator, a deconvolution process may be used to obtain dependable results. 
20 Moreover, the convolving function for a wide aperture depends only on the geometry 
of the collimator and may be expressed as a set of linear equations that can be readily 
solved. 

Thus, in accordance with Ae present invention, a wide-bore collimator probe, 

having a position tracking system, operating at very short time intervals, and 
25 coimected to a data processor for wide-aperture collimation - deconvolution 

algorithms, can be used for the construction of an image of radiation density in one-, 

two-, or three-dimensions, with high degree of accuracy. 

It will be appreciated that similar analyses may be employed for wide aperture 

probes having non-circular collimators. For example, the wide-aperture probe may 
30 have a rectangular collimator of a single cell or of a coarse grid of cells. Additionally, 

the collimator may be a wide-angle collimator. A corresponding deconvolution 

fimction may be used by the data processor. 
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Reference is now made to Figures 27A and 27B, which illustrate a radiation 
sensor 600, preferably generally shaped as a tube collimator Radiation quanta 602 are 
registered by the radiation sensor 600, as described hereinabove, thereby providing the 
average number of quanta per unit time. The radiation sensor 600 may be moved 

5 around a volume of interest 604, preferably, with six-dimensional freedom of motion, 
and six dimensional monitoring. However, it will be appreciated that a more limited 
motion and corresponding monitoring, for example, only in the X and p directions is 
possible. The position of the sensor 600 and it<; direction (as well as the position of the 
investigated volume 604) are assumed to be known at any given moment (Figure 27 A). 

10 The tube collimator is preferably provided with a plane circular detector 606 of 

radiation quanta. The quanta detector 606 is preferably disposed on a rear end 608 of 
the tube and radiation quanta can reach the detector 606 only through an open front 
end 610 of the tube (Figure 27B) 

Reference is now made to Figure 27C, which illustrates a system of coordinates 

15 (x, y, z) with the origin O in the center of the radiation s^or 600, die (x, y) plane 
being the plane of the detector, and the z axis being in the center of the collimator 
tube. The geometry paramet^s of the collimator tube - height h and radius p ~ are 
known. 

From the rotational symmetry of the tube, it is clear that having a radiation 
20 source Q = Q{x, y, z) with the total intensity I uniformly radiating in all directions, the 
portion of the intensity registered by the quanta detector 606 of the radiation sensor 
600 is determined only by the distance r from Q to the axis of the collimator (axis z) 
and the distance z from Q to the (x, y) plane. In other words, there is a function Q(r, z), 
which is defined only by the collimator parameters p and h (corresponding expressions 
25 from p, h, r and z may be easily written in explicit form), such that the intensity of the 
radiation spot Q = Q(x, y , z) == Q(r, z) registered by detector 606 is proportional to 
□{r, z) and to the total intensity I of the radiation spot 

Reference is now made to Figure 27D. It follows from the foregoing 
discussion, that if instead of one radiation spot there is a radiation distribution I(Q) - 1( 
30 Q(r, z) ) in a volume V, then the radiation intensity, registered by radiation sensor 600, 
is proportional (with some constant not depending on the radiation distribution and the 
sensor position) to the following integral: 
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l(Q(r,z))<t>(r,z)dQ (1) 
V 

An algorithm for estimating the intensity distribution I(Q) from the values 
obtained in the measurement scheme of Equation (I) is now discussed. For the sake of 
simplicity, the first case is discussed with reference to Figure 27E for a 
5 two-dimensional problem, wherein intensities I(Q) are distributed in some 
2-dimensional plane. The 3D problem is a direct generalization of the corresponding 
2D problem, as is explained hereinbelow. 

As seen in Figure 27£, the radiation sources are distributed in a rectangular 
region V in a plane. Two systems of coordinates are considered. The 6r$t one is the 
10 sensor coordinate system (x, y» z) corresponding to fhe sensor 600. The second one is 
the radiation source coordinate system (u, v, w) correspomling to the radiation sources 
plane (u, v). 

It is assumed that at each discrete time inciemrat, the position of the origin of 
(x, y, z) system and the direction of the z-axis unit vector in (u, v, w) coordinates are 
15 known. In other words, the position and direction of the moving sensor in the (u^ v, w) 
coordinate system is known, and the (u, v, w) coordinate system is assumed to be 
motionless. 

The radiation sources are considered to be distributed in acconlance with the 
distribution function I(Q) ) in some bounded, given rectangle V on the plane (u, v). 
20 1(0) = I(u, v) is the unknown and sought-for radiation (or radiation intensity) 
distribution function defined in V. 

To regularize the problem of estimation of the radiation distribution function 
I(Q), the fimction I(Q) will be considered to be given li^om some finite dimensional 
space H of functions defined in V. In other words, tiie function I(Q) itself will not be 
25 estimated but rather some finite dimensional approximation of the distribution I(Q). 

The simplest approach to finite dim^sional approximation is to subdivide tlie 
rectangle V into sets of equal rectangular cells and consider the space H of 
step-functions corresponding to this subdivision (i.e., the space of functions that are 
constant in the cells of subdivision), as shown in Figure 27F. 
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If the subdivision of rectangle V into small rectangles is sufficiently fine, then 
this stcp-fimction 2q;)pn)ximation is gpod enou^ for th estimation of radiation 
distribution I(Q). 

Let each side of rectangle V be divided into n equal parts (Figure 27F). Then m 
5 = n^ is the dim^ion of the space H of step-functions on the corresponding 
subdivision. 

The space H is naturally isomorphic to the m-dimensional space of n x n 
matrixes (with its natural scalar product <•,•>). 

LetI = (Iij)i,j - be the unknown element of H that it is desired to 
10 estimate. Suppose that elonent I is measured on K functionals {<I>k } ie==i.x of the form 
of the integral (1): 

<I,<I>k> = 2ij=,..„Iij (2) 

15 where <^ = {O ij i.„ji , k = 1, K (after approximation of function 

J(Q), by the corresponding step-fiinction, the integral (1) is transformed to the sum 
(2)), 

Functionals Ojt, k = 1, K, correspond to K discrete positions of the sensor 
(Figure 27E). Knowing the explicit expressions for functions □(r, z) from (1) and 
20 knowing for each time moment k, the position of the sensor relative to inspection 
region V, one can calculate all matrixes Ok = (^ i j ) i. j - u.^ , k = I , . . K. 
Accordingly the following scheme of measurements are obtained: 
Vk = <I,<I>k> + Ck, k=l,...,K. (3) 
Here \\f k are results of measurements of the unknown element I of the space 
25 H, and 8 k are random errors (c k - independent random variables, Ee k = 0 , k = 1, 
K). 

Let M: H -> H the operator in the space H of the form: 

M -Swx Ok®*k- (4) 
Then the best non-biased linear estimate I of the element I is given by the 
30 formula: 

1 = 1^*4^, <5) 
where M"^ H H the inverse operator to the operator M of the form (4), and 
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= 2k=!..K Vk*k , (6) 
(where k are the results of measurements of the form (3) ). 
One problem of using estimates (5) (besides computational probl^xis if the 
dimension m of the space H is very large) is that flie operator M: H -> H of the form 
5 (3) is "bad invertable". In other words, the estimation probim is "ill-posed". It means 
that having a noise e ^ in the measurements scheme (3), even if the noise is small, may 
sometimes result in a very large estimation error dist(I , I). 

This means that the estimation problem requires additional regulaiization. This 
is a general problem of solving a large set of linear equations. There are several 
10 methods for solving such equations. Below is described one of the known methods for 
solving such equations but numerous other methods are also possible, theses include 
gradient decent methods such as in (http://www-visl.tcchnion.ac.i1/l 999/99-Q3/www/) 
and other methods that are generally known in the art. Further* it is possible to 
improve the image reconstruction by taking into account the correlation between 
15 measurements as they are done with substantial overlap. Also, in the following 
description, a regular step function is assumed for the representation of the pixels or 
voxels, other basis may be used such as wavelet basis, gaussian basis, etc., which may 
be better suited for some applications. 

To obtain regularized estimate I r instead of the estimate t the eigenvector 
20 decomposition of the operator M may be used: 

Let 9i, (p2, . . . , (pm be eigenvectors of operator M: H -> H corresponding to 
eigenvalues A.] > A.2> ... > Xm^O. 

Let R be some natural number, 1 < R < m ( R is the 'Yegularization parameter" 
). Let H^^^ be the subspace of the space H spanned by the first R eigenvectors • . , 
25 <Pr. 

H^'^^ = sp{(pk . (7) 

Let P^^^ : H -> H^*^^ be the orthogonal projection on subspace H^^ . 

The regularized estimate 1 r may be obtained as follows: 

Let Oi,^^^ -P^^Ok, k=l,...,K. 

30 ^'^'^^ = Iic=i..jc (8) 

^(R); ij(R} ^ jj(R) ^jjg operator of the form: 

M^^^ =Xjc^, . ,c .<l>k^^^®a)k<^> (9) 



10/07 03 THU 13:37 FAX 972 3 6127575 



G. E. EHRLICH 



Q04S 



43 

(operator M^^ is the restriction of the operator M of the foim (4) to the 
subspace H^^ of the form (7) ), 

theQi^>-(M^>)"^ (10) 

When the regularization parameter R is properly chosen (so that the eigenvalue 
5 Xr is not too small), fiien the estimate (10) becomes stable. 

There are several possible approaches to choosing the parameter R. One 
approach is to leave R as a ''program parameter'' and to obtain the reasonable value ''in 
experiment**. Another appn>ach is to choose some "optimal" value. This is possible if 
the covariation operators of the random noise 8 k in (3) are known, and information 
10 about the element I of the space H is known a priori. 

The subdivision of rectangle V into a large number of equal rectangles has the 
disadvantage of making the dimension of the space H too big (especially in the 3D 
case). If each side of rectangle V is subdivided into n equal parts, then the dimension 
of the space H will be and the dimension of the matrices used in solving the 
15 corresponding estimation equations would be n^ x = n'* (in the 3D case, n^ x = 
n^. It is clear that for large n, this situation may cause serious memory and 
computation time. 

In accordance with a preferred embodiment of the present invention, an 
irregular subdivision of Ihe rectangle V is used. This irregular approach may 
20 significantly decrease the dimension of the problem and faciUtate computer 
calculations. 

More specifically, a drawback of the regular subdivision of the investigated 
region V, discussed hereinabove, is that a lot of cells that actually have no signal may 
be taken into account (Figure 27F). It would be much better to have small cells only in 
25 regions with high signal and have big cells in regions with low signal. 

Reference is now made to Figure 27G, which illustrates an advantageous 
irregular cell subdivision of rectangle V, in accordance with a preferred embodiment 
of the present invention. 

In a first stage, regular subdivisions arc made in "large" cells, and 
30 measurements and estimations are made as described hereinabove. In this manner, the 
intensity distribution is estimated in the large cells. 
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In a second stage, the large cells> which have an intensity larger than some 
threshold, are subdivided into 4 equal subcells (or 8 subcells in the 3D case). A 
suitable threshold may be obtained by taking the average intensity (of all large cells) 
minus two (or three) sigmas (standard deviation^ for example. Measurements and 
S estimations are made in these subdivisions as described h^einabove. 

The act of subdivision and subsequent measurements and estimations are 
continued until a desired accuracy is reached at some smaller level of subdivision, 
typically defined by the computational and memory capabilities of the compute- being 
used 

10 The 3D problem may be treated in the same way as the 2D case, fhe only 

difference being that instead of rectangle V, there is a parallelepiped V (Figure 27D). 
Accordingly, the cells in each subdivision are also parallelepipeds. 

In accordance with the' present invention, the wide-aperture collimatioa - 
deconvolution algorithms described hereinabove may be used for a variety of imaging 

15 systems. For example, the algorithms may be used with single radiation detector 
probe, an array of radiation detector probes, large gamma cameras of various designs, 
sudi as multi-head cameras, genera] purpose cameras, and automatic white balance 
(AWB) scanners. The algorithms are suitable for SPECT and planar imaging, and may 
be used for all types of isotopes of with any type of photon energy. 

20 From the foregoing discussion, the skilled artisan will appreciate that the 

algorithms described hereinabove may be used to predict the location of the radiation 
source and the uncertainty region (based on the system measurement errors) in the 
vicinity of the radiation source. The algorithms also guide the user to perfom 
additional measurements to minimize the uncertainty region according to the 

25 requirements of the system operator. 

The algorithms thus comprise a feedback system that employs analysis to 
determine the bounds of an uncertainty region about the radiation source, and which 
guides medical persormel to conduct additional scans in these uncertainty regions to 
improve accuracy, reduce error, and hence minimize the bounds of the uncertainty 

30 regions. 
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Continuous sampling with rddialion probe 202 may provide localization of a 
tumor and a physiological radiation activity map of the tumor regiuiu Higher safety 
and accuracy are gained by a greater number of scans. 

Other deconvolution methods are known and are often used in image 
5 processing procedures. Examples of such deconvolution methods are described in US 
Patent 6,166,853 to Sapia et al., the disclosure of which is incorporated herein by 
reference. (However, it is appreciated that these are just examples and the present 
invention is not limited to the deconvolution methods mentioned in US Patent 
6,166,853) 

10 Reference is now made to Figure 22 which illustrates a flow chart of a 

radiation map reconstruction algorithm, by another system of wide-aperture 
coUimation - deconvolution algorithms, in accordance with another preferred 
embodiment of the present invention. 

In typical image acquisition, Hght (or other electromagnetic wave energy) 

1 5 passes through a jGnite aperture to an image plane. The acquired image is a result of a 
convolution of the source object's ligiht with the aperture of the imaging syston. A 
system transf^-function may be generally obtained directly by taking the Fourier 
transform of the aperture. As is known in the art, the blurring effects due to 
convolution generally exist in two-dimensions only, i.e., the x-y planes. A 

20 point-spread-fiinction (PSF) is an expression used to describe the convolutional 
blurring in two-dimensions. The PSF physically results from imaging a point source. 
The Fourier transform of the PSF is the system transfer-ftmction, which is obtained by 
convolving the system transfer-function Avitli a Dirac-deita function. A point source is 
the physical equivalent of a Dirac-delta function, and, in the frequency domain, the 

25 Dirac-delta function is a unity operator across the spectrum. Therefore, the Fourier 
transform of the PSF should be the Fourier transform of the aperture. However, the 
PSF contains noise and blurring due to other effects such as aberrations. 

The PSF contribution to the overall blurriness may be diminished or eliminated 
by deconvolution. 

30 In the case of the present invention, the transfer function of the radiation 

detector may be determined by taking the Fourier trarksform of the aperture of the 
detector, and taking into account the noise and blurring due to other effects such as 
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abeirations (step SOO). An example of a transfer function may be a normal distribution. 
Using known mathematical techniques, the deconvohition of the transfer function may 
be detennined (step 502). 

The count readings of each spatial location of the detector constitute the sum of 

5 radiation counts from all the voxels (or pixels in the case of two-dimensional maps» the 
term "voxel" being used herein to include both pixels and voxels) within the detector's 
field of view. At least one voxel, or preferably each such voxel, may be assigned a 
count value based on the deconvolution of the unique transfer function of the radiation 
detector in use (step 504). An additional mathematical procedure may treat the various 

10 values that each voxel receives due to the multiple readings from viewpoints of 
different detectors (step 506). This treatment may constitute for example a simple 
algebraic aveiage, minimum value or reciprocal of averaged reciprocals in order to 
produce a single value of readings in each voxel. The deconvolution is then used to 
reconstruct the voxels of the radiation map with diminished or no bluniness (step 

15 508). 

The algorithms described herein are applicable not only to the analysis of 
readings obtained using a directional radioactivity detector, rather they also apply for 
spatially sensitive (pixelated) radioactivity detectors, hi this case, the readings of each 
pixel are algorithmically treated as described herein like for a directional radioactivity 

20 detector. The motivation bdund using a spatially sensitive detector is to save on 
measurement time by receiving readings from a multitude of directions in parallel. 
This, in essence, creates a number of overlapping low resolution images which can 
then be processed to form a high resolution image. In addition, the spatially sensitive 
detector can be scanned to improve even further the resolution using the algorithms 

25 described hereinabove. 

Thus, the same algorithms that apply for a directional detector apply for the 
spatially sensitive detector, only now instead of one radiation reading at each position, 
a large set of desecrate positions are processed in parallel. Each pixel can be seen as a 
separate detector with an angle of acceptance dictated by the geometry of a segmented 
30 collimator employed th^eby. Each of the pixels occupies a different position in space 
and hence can be seen as a new position of a single directional probe by the algorithm 
described herein. It is also possible, like with the directional detector, to scan the 
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whole set of pixels by scanning the spatially sensitive detector and to acquire a new 
set of data points from the new position. Once obtaining a low resolution image finom 
each of the pixels of the spatially sensitive detector, a super resolution algorithm can 
be employed to generate an image of higher resolution. Suitable super resolution 
5 algorithms are described in, for example, J. Acoust Soc. Am., Vol. 77, No. 2, 
February 198S Pages 567 - S72; Yokota and Sato, IEEE Trans. Acoust. Speech Signal 
Process. (April 19S4); Yokota and Sato, Acoustical Imag^g (Plenum, New York, 
1982, Vol. 12; H. Shekarforoush and R. Chellappa, "Data-Ehiven Multi-channel 
Super-resolution with Application to Video Sequences'*, Joumai of Optical Society of 

10 America-A, vol. 16, no. 3, pp. 481-492, 1999; H. Shekarforoush, J. Zenibia and M. 
Berthod, '^Extension of Phase Correlation to Sub-pixel Registration", lEfiE Trans. 
Image Processmg, to appear, P. Cheeseman, B. Kanefsky, R. Knift, J. Stutz, and R. 
Hanson, "Super-Resolved Surface Reconstruction From Multiple Images," NASA 
Technical Report FIA.94-12, Dec«nber, 1994; A.M. Tekalp, M.K. Ozkan, and M.L 

15 Sezan, "High-Resolution Image Reconstruction for Lowcr-Resohition Image 
Sequences and Space-Varying Image Restoration," IEEE Xntemational Conference on 
Acoustics, Speech, and Signal Processing (San Francisco, CA), pp. HI- 169- 172, 
March 23-26, 1992, http://www-visl.technion.ac.il/1999/99-03/www/. w hich are 
incoiporated herein by reference. 

20 Referring fimher to the drawings, Figures 28A - 28F schematically illustrate a 

radioactive emission probe 700, preferably, handheld, for free-hand scanning, in 
accordance with preferred embodiments of the presoit invention. 

As seen in Figure 28A, handheld probe 700 is a lightweight nuclear imaging 
camera, adapted for fice-hand movement and having a positioning device 714 

25 incoiporated thereto. 

Handheld probe 700 includes a radiation detector 706, mounted onto a housing 
702, wherein the various electronic compouCTts (not shown) are contained. Housing 
702 is preferably formed of a rigid, lightweight plastic, a composite, or the like, and 
includes a handle 704, for easy maneuvering of probe 700. 
30 A control unit 710, which may be mounted on housing 702, may include basic 

control knobs, such as "stop", "start", and "pause." Additionally or alternatively, 
control unit 710 may be a computer unit, such as a microcomputer or the like, having 
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processing and memory units, contained within housing unit 702. Additionally, 
control unit 710 may include a display screen 718, mounted on housing unit 702, for 
displaying bfonnation such as gamma counts, gamma energies, device position, and 
the like. Display screen 718 may be interactive, for control and display. Control unit 
5 710 may further include a data unit 720, for receiving a diskette, a minidisk, or the 
like. Preferably, data unit 720 is a read and write unit. An eject button 722 may be 
included with data unit 720. 

A cable 724 may provide signal communication between handheld probe 700 
and a main-data-collection-and-analysis-system 726, such as a PC computer or a 

10 server. Additionally or alternatively, signal communication with 
main-data-coUection-^and-analysis-systc^ 726 may be wireless, for example, by radio 
frequency or infrared waves. Alternatively, hand-held probe 700 may be a stand-alone 
unit, optative with built-in con:q>uter unit 710. Additionally, data may be transferred 
to and from main-data-coUection-and-analysis-system 726 via a diskette or a minidisk, 

15 from data unit 720. 

A preferably rechargeable battery 716 may be attached to housing unit 702. 
Additionally or alternatively, a cable 712 may be used to provide power 
communication between handheld probe 700 and the grid (not shown), or between 
handheld probe 700 and main-data-collection-and-analysis-system 726. 

20 Positioning device 714 is preferably a Navigation sensor, for determination of 

six coordinates X, Y and Z axes and rotational angles p, 8 and (p, as shown in Figure 
28F. Positioning device 714 may be mounted on housing 702, or located within it. 
Preferably, positioning device 714 is a magnetic tracking and location system, known 
as miniBBEU^, commercially available from Ascension Technology Corporation, P.O. 

25 Box 527, Burlington, Vermont 05402 USA 

(http://www.ascension-tech.coro/graphic.htm). Alternatively, any other known, 
preferably 6D tracking and positioning system may be used. Positioning device 714 
maintains signal communication with either one or both of control unit 710 and 
main-data-collection-and-analysis-system 726. 

30 As seen in Figures 28 A and 28B, nuclear imaging detector 706 may be a disk, 

forming a single pixel, having a diameter D of about 10 min, and a length LI of about 
5 mm. It will be appreciated that oth^r dimensions, which may be larger or smaller. 
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are also possible. Preferably, radiation detector 706 is connected to a preamplifier 
(not shown). 

Additionally, as seen in Figures 28A and 28B, handheld probe 700 may 
include a tubular, wide-bore collimator 708, fonncd for example, of lead or of 

5 tungsten. The collimator may have a length L2 of about 10 to 30 mm, and preferably, 
about 15 mm. It will be appreciated that other dimensions, which may be larger or 
smaller, are also possible. 

Alternatively, as seen in Figures 28C and 28D, collimator 708 may be a square 
grid of 4 X 4 ceUs 705, or 16 X 16 cells 705. Collimator grid 708 m^ have sides w 

10 of about 1 0 mm, and length L2 of about 20 - 40 mm. It will be appreciated that other 
dimensions, which maybe larger or smaller, are also possible. 

As seen in Figure 28E, each cell 705 may mclude a plurality of radiation 
detector pixels 703, formed for example, as square pixels of 5 X 5 mm and length LI 
also of bout 5 mm. The reason for the division of nuclear detector 708 into pixels 703 

IS is that while the wide aperture of a coarse collimator increases the counting efficiency, 
radiation detector uniformity is better maintained when the size of the radiation 
detectors is small. Thus, it is desired to combine a coarse collimator grid with a fine 
pixel size. Lead spacers are placed between the pixels. From spatial resolution 
consideration, each cell 70S preferably operates as a single pixel, but each pixel 703 

20 within cell 705, must be calibrated individually, to conrect for material nonuniformity, 
as will be described hereinbelow, in conjunction with Figure 29B. 

Preferably, radiation detector 706 is equipped with a positive contact array, for 
forming contact with each pixel 703, and a common negative contact Preferably, 
each pixel 703 is connected to a preamplifier. 

25 Preferably, radiation detector 706 is a single module anray, for example, of 4 X 

4, or of 16 X 16 pixels, of room temperature CdZnTe, obtained, for example, from 
IMARAD IMAGING SYSTEMS LTD., of Rehovot, ISRAEL, 76124, 
wM^.imarad.com . 

Room temperature solid-state CdZnTe (CZT) is among the more promising 
30 nuclear detectors currently available. It has a better count-rate capability than oth^ 
detectors on the market, and its pixilated structure provides intrinsic spatial resolution. 
Furthermore, because of the direct conversion of the gamma photon to charge-cairiers, 
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energy resolution is enhanced and there is better rejection of scatter events and 
improved contrast. 

In accordance with the present invention, the detector may be optimized, in 
accordance with the teaching of '"Electron lifetime determination in semiconductor 
5 gamma detector arrays " http://urila.tripod.com/hechthtm> "GdTe and CdZnTe 
Crystal Growth and Production of Gamma Radiation Detectors," 
http://members.lripod.com/'-urila/crvstaLhtm, and "Driving Energy Resolution to the 
Noise Limit in Semiconductor Gamma Detector Arrays/' Posto- presented at 
NSS2000 Conference, Lyon France, 15 - 20 October 2000. 

10 http://urila.triDod>com/NSS.htm. all by Uri Lachish, of Guma Science, P.O.Box 2104, 
Rehovot 76120, Israel, urila@intemet-zahav.net, all of whose disclosures are 
incorporated herein by reference. 

Accordingly, radiation detector 706 may be a monolithic CdZnTe crystal, 
doped with a trivalent donor, such as indium. Alternatively, aluminum may be used 

15 as the trivalent donor. When a trivalent dopant, such as indium, replaces a bivalent 
cadmium atom within the crystal lattice, the extra electron falls into a deep trap, 
leaving behind an ionized shallow donor. The addition of more donors shifts the 
Fermi level from below the trapping band to somewhere within it. An optimal donor 
concentration is achieved when nearly all the deep traps become occupied and the 

20 Fermi level shifts to just above the deep trapping band. 

Optimal spectral resolution may be achieved by adjusting the gamma charge 
collection time (i.e., the shape time) with respect to the electron transition time from 
contact to contact. Gamma photons are absorbed at different depth within the detector 
where they generate the electrons. As a result, these electrons travel a different 

25 distance to the counter electrode and therefore produce a different external signal for 
each gamma absorption event. By making the shape time shorter than the electron 
transition time, from contact to contact, these external signals become more or less 
equal leading to a dramatic improvement in resolution. 

Furthermore, for a multi-pixel detector, the electrons move from the point of 

30 photon absoiption towards the positive contact of a specific pixel. The holes, which 
arc far slower, move towards the negative contact, and their signal contribution is 
distributed over a number of pixels. By adjusting the gamma charge collection time 
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(i.e.» the shape time) with respect to the election transition time from contact to 
contact, the detector circuit collects only the electrons' contribution to the signal, and 
the spectral response is not deteriorated by the charge of the holes. 

For an optimal detector, crystal electrical resistively may be, for example, 

5 about 5 X 10 ^ ohm cm. The bias voltage may be, for example, - 200 volts. The 
shape time may be, for example, 0.5 X 10 ~^ sec. It will be appreciated that other 
values, which may be larger or smaller are also possible. 

In accordance with a preferred embodiment of the present invention, control 
unit 710 includes data acquisition and control components, which further include a 

10 veiy small, single-module Detector Carrier Board (DCB). The DCB preferably 
includes a temperature control system, having a temperature sensor, a desired 
temperature setting, and a heat removal system, preferably operative by peltier 
cooling. Preferably, the DCB further includes a High Voltage (HV) power supply, a 
HV setting, a HV and dark current output-sensing, and power regulators, preferably 

IS having a range of about i:2V, a XAIM MBAIS cuncnt driver, and noise filtering 
components. 

Preferably, data acquisition and control components are based on Application 
Specific hitegrated Circuit (ASIC) architecture for the control, calibration, and 
readout. ASIC data acquisition and control includes calibration of ASIC data, W 

20 conversion, amplification and shaping. Data identification includes threshold 
discrimination (i.e., noise rejection) and energy windows discrimination, wherein 
tha*e may be one or several energy windows per radionuclide. Preferably, there are up 
to four energy windows per radionuclide, but it will be appreciated that another 
number is also possible. 

25 Data acquisition relates to counting events per pixel, per energy window, for 

an acquisition period, which may be very short, for example, about 100 msec, in a 
continuous mode. It will be appreciated that larger or smaller counting periods arc 
also possible. 

The events per pixel 703, per energy window, is synchronized with signals 
30 fix>m positioning device 7 14, preferably via a DIO card, either within control unit 710, 
or main-data collection-and-analysis-system 726, or both. Additionally, initialization 
information, including Digital to Analog Converter (DAC) information, may be stored 
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in software of main-data-collection-and-analysis-system 726 or within control system 
710. 

Preferably, an Integrated Drive Electronics (IDE) system provides the 
electronic interface. 

s Referring further to the drawings, Figures 29A - 29B schematically illustrate 

the manner of calibrating handheld probe 700, in accordance with preferred 
embodiments of the present invention. Non-uniformity of the material leads to 
different sensitivies for each pixel. A correction is required to eliminate the 
sensitivity effect of the different pixels. 

10 An energy spectrum for each pixel 703 of radiation detector 706 is obtained in 

a classical fashion. Radiation detector 706 is irradiated by a flood source of a known 
photon energy, for example, 24 KeV, and for each pixel 703, such as 703(1), 703(2), 
703(3)... 703(n), a summation of counts as a function of voltage is obtained, during a 
finite time interval. The location of the energy peak, corresponding to the peak 

15 photon energy, for example, 24 KeV is bounded, by an upper level (UL) and a lower 
level (LL) for that pixel 703(n), forming an energy window for the specific pixel. 
Only photons which fall within the bounded energy window for each pixel are 
counted, to exclude Compton scattmng and pair production. A multi-energy source 
may be used, and several energy windows may be bounded for each pixel. 

20 After the energy windows, representing specific energy peaks, are detmnined 

for each pixel, they are stored in a memory unit. 

During acquisition, each collimator cell 705(n) is operative as a single pixel. 
Thus, the sensitivity correction factor has to be at a collimator-cell level, rather than at 
a pixel- level. 

25 To calculate the sensitivity correction factor, an avCTage counts per cell 705, is 

obtained as: 



tt ^r.^ sum of counts in all pixels 703 of detector 706 
average counts per cell 705 = — — 

number of collimator cells 705 in detector 706 
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correction factor for cell 705(ii) = average counts per cell 705 

counts in a cell 705(n) 

During aquistion, the sum of counts within the bounded windows of all pixels 
in cell 705(n) is multiplied by the the correction factor for that cell 705(n) to obtain a 
sensitivity corrected count rate. 

5 Refenring further to the drawings. Figure 30 schematically illustrates the 

manner of synchronizing event readings by radiation detector 706 and position 
readings by positioning device 714, in accordance with preferred embodiments of the 
present invention. Preferably, event readings, for each cell 705 is tabulated with the 
position readings for each time interval, such as t2> t3, after correction for 

10 sensitivity. 

Referring further to the drawings, Figures 31A - 31 C describe a spatial 
resolution bar-phantom test of probe 700, in accordance with a preferred embodiment 
of the present invention. As in Figures 17 and 18, hereinabove, the aim of the bar 
phantom test is to evaluate image quality by measuring image resolution of a bar 
15 phantom. Spatial resolution is expressed in terms of the finest bar pattern, visible on 
nuclear images. 

As seen in Figure 31 A, a bar phantom 750 was constructed of lead strips 752, 
encased in a plastic holder (not shown). The lead strips were of equal width and 
spacing. Thus, a 2- mm bar pattern consisted of 2-mm wide strips 752, separated by 
20 2-mm spaces 754, along an x axis. The thickness of lead strips 752 was, sufficient to 
ensure that bars 752 are virtually opaque to the y rays being imaged. 

A uniform radiation field, formed of a sheet source (not shown), was placed 
directly behind the bar phantom 750, and bar phantom 750 was scarmed with probe 
700. 

25 Probe 700 of the present embodiment comprised an eV Pen-type CdZnTe soUd 

state detector with collimator diameter D, of 5 mm, and collimator length, L2 of 30 
nun. The CdZnTe detector was attached to a charge sensitive preamp probe and 
connected to a Digital and analog V-target alpha system PC cards. The probe energy 
window calibration was performed at the beginning of the test by using EG&G Ortec 

30 Maestro 2K Tramp ISA Spectra Processing connected to the analog and the digital 



10/07 03 THU 13:42 FAX 972 3 8127575 



G. E. EHRLICU 



@056 



54 

V-target alpha system cards. The probe position infonnatiou was provided by 
MicroScribe - 6DLX arm manufactured by Immersion Coiporation. Acquisition time 
was 6 minutes. 

Figure 3 IB illustrates a bar phantom image 760, produced by probe 700, after 
5 6 minutes of acquisition time, in accordance with the present invmdon. Image 760 
shows bar images 756, and space images 758, of the 2-mm bar pattern, clearly 
demonstrating a hig|h degree of spatial resolution. 

Figure 31C provides a similar observation, based counts as a function of 
length axis, X. Tlie well-defined, narrow sinusoidal peaks of counts as a function of 
10 length axis demonstrate a hig}i degree of ^atial resolution for the 2-mm bar pattern. 

Referring further to the drawings, Figures 32A - 32D describe a spatial 
resolution bar-phantom test of a prior-art probe, a High Resolution Gamma Camera of 
a competitor. Although for the prior art probe, acquisition time was 3 hours, a bar 
pattern is barely observed for a 3-mm bar pattern of Figure 32C, and clearly observed 
15 only for a 3.5 mm bar pattern of Figure 32D. Table 1 provides a comparison with a 
prior art technology. 



Table 1 



Imaging 

Technology 


Detector - 
collimator 

arrangement 


Acquisition Time 


Maximum Spatial 
Resolution 


a 

High-resolution 
Gamma-Camera 
of a 
competitor's 


multi-pixel 
detector, 
tube collimators 
1.9ramindia. 
40 mm in 
length 


180 min. 


Worse then 
3 mm 

• 


Probe 700 
of Figure 28 A 


single-pixel 

detector 
10 mm in dia. 
tube collimator 
25 mm in length 


* 

6 min. 


2 mm 


Summary 




V-target 
probe is 30 time 
faster 


V-target's 
resolution is 2 time 
better 



As seen in Table I, there is a clear advantage to the probe of the present 



invention. 

20 Good energy resolution of the detector is important for precise identifica.tion 

and separation of y rays, for radionuclide identification and for scatter rejection. A 
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major advantage of CdZnTe detector is better energy resolution, for example than that 
oftbeNai(Tl) detector. 

Referring fiuther to the drawings. Figures 33A -33B illustrate the energy 
resolution of a single pixel 703 of probe 700, . in accordance with the present 
5 invention. 

Figure 33A compares a spectrum of a single pixel of a CdZnTe detector with 
that obtained by a single pixel of an NaI(Tl) detector, at ~ 120 KeV, showing the 
difference in enCTgy resolution between the two detectors. Clearly, the a CdZnTe 
detector provides a sharper, more well-defined energy peak, which is about 1/3 higher 
10 than the Nai(n) peak, and having a FWHM that is 2/3 that of the Nai(Tl) peak. 

Figure 33B provides a spectrum ofaCdZnTe detector, showing good energy 
resolution at the low energy level of about 24 KeV. 

Referring further to the drawings. Figures 34A - 34C schematically illustrate 
endoscopic radioactive emission probes 800, adapted to be inserted into the body, on a 
15 shaft or a cathet^ 713, via a trucar valve 802 in a tissue 810. 

As seen in Figure 34A, endoscopic radioactive emission probe 800 may 
include, for example, radiation detector 706, for example, of a single pixel, wide-bore 
collimator 708, position tracking system 714, and an extracorporeal control unit 804. 
Preferably, the diameter of radiation detector 706, Ddettetor> is less than 12 mm. In 
20 accordance with another embodiment of the present invention, radioactive emission 
probe SOO may be used in body lumens, for example, for scanning the stomach or the 
uterus. 

As seen in Figure 34B, mdoscopic radioactive emission probe 800 may 
include several radiation detectors 706, for example, of a single pixel each, and having 
25 a dedicated collimator 708, each. Unhke the collimator grid structure of Figures 28C 
and 28D, the collimators of Figure 34B are not parallel to each other, and each may 
point m a different direction. Alternatively, each radiation detector 706 may include a 
plurality of pixels. 

As seen in Figure 34C, endoscopic radioactive emission probe 800 may 
30 include several radiation detectors 706, which may include wide-angle collimators. 

Referring fiirther to the drawings, Figure 35 schematically illustrates a manner 
of calculating a distance, d, between radioactive emission probe, such as probe 700 
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(Figures 28 A - 28F) and a radiation source 766, based on the attenuation of photons of 
different energies, which are emitted from the same source. Probe 700 is on an 
extracorporeal side 762 and radiation source 766 is on an Ihtracorporeal side 764 of a 
tissue 768. 

The transmitted intensity of a monochromatic gamma-ray beam of photon 
energy E and intensity Jq crossing an absorbing layer, such as a tissue layer, of a depth 
d, is a function of the photons energy £, the layer thickness d, and the total linear 
attenuation coefficient as a function of energy, ^totai> of the material constituting the 
layer. The transmitted intensity is given by: 

f(E) = lofE) expf-ii^) d) (Eq-I) 



The linear attenuation coefficient is expressed in cm'^ and is obtained by 
multiplying the cross section (in cm /gr) with the absorbing medium density. The total 
1 5 linear attenuation coefficient accounts for absorption due to all the relevant intmction 
mechanisms: photoelectric effect, Compton scattering, and pair production. 

The layer thickness d, which is the depth of the radiation source, can be 
measured using a radionuclide that emits two or more photon energies, since 
attenuation is a function of the photon energy. 
20 For photons of different energies, originating fiom a same source, and having 

the same half-life, we get: 

m) I,(E,} e->'''^>' 



^^A^ [HE,)/ HE,)] 



(Eq-4) 
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5 TThe source depth, d, can be calculated as: 



The linear attenuation coefficients for specific energies are calculated from the 
10 photon mass attenuation and energy absorption coefficients table. 

The ratio Io(E]) / loCEi) for a known isotope or isotopes is calculated from an 
isotope table and the ratio I(E|) / I(E2) is measured at the inspected object edges. 

Additionally, X-ray information, for example, mammogr£q)hy may be used, for 
more material linear attmuation coefficient, since X-ray imaging provides information 
1 S about the material density. 

For example, we may consider a lesion inside a soft tissue, such as a breast, at 
a distance d that holds an isotope of ^^^/, which emits two main photons, as follows: 
El = 27kcV with 86.5% ;and 
E2=159keVwi(h 83.4%. 
20 The soft tissue density is 0.8g/cm\ 

The linear attenuation coefficients for Ej and E2 are calculated from J.H 
Hubbell, "Photon Mass Attenuation and Energy - absorption Coefficients from IkeV 
to 20MeV", Int. J. Appl. Radiat. Isat. Vol. 33. pp. 1269 to 1290, 1982, as follows: 
H(E,) = 0,474 cm/g -> ^(Ei) = 0.474 cm/g • 0.8 g/cm^ = 038 1/cm 
25 ^(£2) = 0.1489 cm/g -> |i(E2) = 0.1489 cm/g • 0.8 g/cm^ = 0.119 1/cm 

The measured ratio between I{Ei) to I(E2) on the attenuating soil tissue edge is 

0.6153. 

Using Eq-7, one gels that the lesion depth is: d 2 cm. 
Alternatively, if the measured ration between I(Ei) to I(E2) is 0,473 the 
30 extracted lesion depth is: d = 3 cm. 
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It will be appreciated that when two or more isotopes used for obtaining the 
different photons, corrections for their respective half-lives may be required. 

It will be ^reciated that by calculating the depth of a radiation source at each 
position^ one may obtain a radiation source depth map» in effect^ a three-dimensional 
5 iamge of the radiation source. This information may be superimposed on the radiation 
source image, as produced by other methods of the present invention, as an 
independent check of the other methods. For example, a radiation source image 
produced by the wide-aperture collimation ~ deconvolution algorithms may be 
superimposed on a radiation source image produced by depth calculations of based on 
10 the attenuation of photons of different energies. 

EXPERIMENTAL RESULTS 

In a series of clinical experiments, some of the basic concepts of the invention 

have been tested on patients who were pre-injected with a suitable radiopharmaceutical 

15 for their particular pathology. Two-dimensional color-coded maps have been 

constructed based on a scan of a pre-detennined lesion area by a hand-held radiation 

detector with a magnetic position-tracking system. The resulting maps, which 

represented the radiation count level, were compared to images of conventional 

18 

gamma camera. The list of radiopharmaceuticals tested includes FDG, 
20 ^^^Tc-MDP, ^^^c sodium pertechnetate, ^^^Tc erythrocytes. Similar 

radiolabeled patterns were observed in the images produced by the system of the 

invention and in the images produced by a conventional gamma camera in the 

following pathologies: 

Figures 23A and 23B illustrate radiolabeled patterns observed in images 
25 produced by the system of the invention and by a conventional gamma camera, 

respectively, of an autonomous adenoma of a thyroid of a 58 year-old male. 

Figures 24A and 24B illustrate radiolabeled patterns observed in images 

produced by the system of the invention and by a conventional gamma camera, 

respectively, of suspected Paget*s disease of a humerus in an 89 year-old female. 
30 Figures 25A and 2SB illustrate radiolabeled patterns observed in images 

produced by the system of the invention and by a conventional gamma camera, 

respect! vely, of chronic osteomyelitis in a 19 year-old female. 
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Figures 26A and 26B illustrate radiolabeled patterns observed in images 
produced by the system of the invention and by a conventional gamma camera, 
respectively, of skeletal metastasis &om inedulloblastoma in an 18 year-old male. 

Figure 36 illustrates a control volume 701, and a three-dimensional image 703 

5 of the radioactivity emitting source, produced by a free-band scaiming of a cancerous 
prostate gland, ex vivo, with band-held probe 700 of Figure 28A. The diameter of the 
detector was 10 nmi, and the length of the collimator was 25 mm. The wide-aperture 
collimation - deconvoluting algorithm was used, for image reconstruction. A scale 
705 illustrates brightness as a function of count rate density. The brighter the area, the 

10 higher the count density. The Figure shows the capability of hand-held probe 700 to 
detect cancerous tissue. The By comparison, the competitor's High-resolution, 
Gamma-Camera failed to detect any cancerous growth. 

The following provides a list of known procedures which can take advantage of 
the system and method of the present invention: 

15 In cancer diagnosis the syste m and method of the present invention can find 

uses for screening for cancer and (or) directing invasive diagnosis (biopsies) either 
from outside the body or by way of endoscopic approach. Examples include, but are 
not limited to, lung cancer biopsy, breast cancer biopsy, prostate canc^ biopsy, 
c^vical cancer biopsy, liver cancer bippsy, lymph node cancer biopsy, thyroid cancer 

20 biopsy, brain cancer biopsy, bone cancer biopsy, colon cancer biopsy, gastio intestine 
cancer endoscopy and biopsy, endoscopic screening for vaginal cancer, endoscopic 
screening for prostate cancer (by way of the rectum), endoscopic screening for ovarian 
cancer, (by way of the vagina), endoscopic screening for cervical cancer (by way of the 
vagina), endoscopic screening for bladder cancer (by way of the urinary track), 

25 endoscopic screening for bile cancer (by way of the gastrointestinal track), screening 
for lung cancer, screening for breast cancer, screening for melanoma, screening for 
brain cancer, screening for lymph cancer, screening for kidney cancer, screening for 
gastro intestinal cancer (from the outside). 

In the special case of MRI, the radiation detector can be combined and 

30 packaged together witihi a small RF coil for the transmission and reception or reception 
only of the MRI signals in a rectal probe configuration for prostate diagnosis and 
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treatment or aiiy other close confinement position such as the vagina, airways, the uper 
portion of the gastrointestinal traclc, etc) 

Procedures known as directing localized treatment of cancer can also benefit 
fiom the system and method of the present invention. Examples include, but are not 

5 limited to, intra tumoral chemothers^y, intra tumoral brachytherapy, intra tumoral 
cryogenic ablation, intra tumoral radio fi^uency ablation, intra tumoral ultrasound 
ablation, and intra tumoral laser ablation, in cases of» for example, lung cancer, breast 
cancer, prostate cancer, cervical cancer, liver cancer, lymph cancer, thyroid cancer, 
brain cancer, bone cancer, colon cancer (by way of endoscopy through the rectum), 

10 gastric cancer (by way of endoscopy through die thorax), thoracic cancer, small 
intestine cancer (by way of endoscopy through the rectum or, by way of endoscopy 
through the thorax), bladder cancer, kidney cancer, vaginal cancer and ovarian cancer. 

hi interventional cardiology the following procedures can take advantage of the 
present invention wherein the method and system can be used to assess tissue 

IS perfusion, tissue viability and blood flow intra operatively during PTCA procedure 
(balloon alone or in conjunction with the placement of a stent), in cases of cardiogenic 
shock to asses damage to the heart, following myocardial infarct to asses damage to 
the heart, in assessing heart failure condition tissue in terms of tissue viability and 
tissue perfusion, in intra vascular tissue viability and perfusion assessment prior to 

20 CABG operation. 

The radioactivity detector can be mounted on a catheter that is entered through 
the blood vessels to the heart to evaluate ischemia fhrni within the heart in order to 
guide ablation probes or another type of treatment to the appropriate location within 
the heart. Another apphcation which may benefit from the present invention is the 

25 localization of blood clots. For example, a radioactivity detector as described herein 
can be used to asses and differentiate between new clots and old clots. Thus, for 
example, the radioactivity detector can be placed on a very small caliber wire such as a 
guide wire that is used during PTCA in order to image intrabloodvessel clots, 
hitrabloodvessel clots can be searched for in the aortic arc as clots therein are 

30 responsible for about 75 % of stroke cases. 

Using the method and system of the present invention to assess tissue 
perfusion, tissue viability and blood flow intra operatively can also be employed in the 
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following: during CABG operation to asses tissue viability, to mark in&rct areas, 
during CABG operations to asses the success of the re vascularization. 

The present invention has many other applications in the direction of 
therapeutics, such as, but not limited to, implanting brachytherapy seeds, ultrasound 
5 microwave radio-fiequency cryotherapy and localized radiation ablations. 

It will be appreciated that many other procedures may also take advantage of 
the present invention. 

It is appreciated that certain features of the invention, which are, for clarity, 
described in the context of separate embodiments, may also be provided in 
10 combination in a single embodiment. Conversely, various features of the invention, 
which are, for brevity, described in the context of a single embodiment, may also be 
provided separately or in any suitable sut>combihation. 

Although the invention has been described in conjunction with specific 
embodiments thereof, it is evident that many alternatives, modifications and variations 
1 S will be apparent to those skilled in the art. Accordingly, it is intended to mibrace all 
such alternatives, modifications and variations that fall within the spirit and broad 
scope of the appended claims. All publications in printed or electronic form, patents 
and patent appUcations mentioned in this specification are herein incorporated in their 
entirety by reference into the specification, to the same extent as if each individual 
20 publication, patent or patent application was specifically and individually indicated to 
be incorporated herein by reference. In addition, citation or identification of any 
reference in this application shall not be construed as an admission that such reference 
is available as prior art to the present invention. 



